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1.  INTRODUCTION 


This  workshop  on  "Dynamical  Instabilities  in  Homogeneously  Broadened 
Lasers"  represents  the  ninth  of  a  series  of  intensive  academic/  government 
interactions  in  the  field  of  advanced  electro-optics,  as  part  of  the  Army 
sponsored  University  Research  Initiative.  By  documenting  the  associated 
technology  status  and  dialogue  it  is  hoped  that  this  baseline  will  serve  all 
interested  parties  towards  providing  a  solution  to  high  priority  Army 
requirements.  Responsible  for  program  and  program  execution  are 
Dr.  Nicholas  George,  University  of  Rochester  (ARO-URI)  and  Dr.  Rudy  Buser, 
NVEOC. 


2.  SUMMARY  AND  FOLLOW-UP  ACTIONS 


Dr.  Pinto  opened  the  meeting  giving  an  overview  of  current  projects  and 
problems  of  interest  in  the  Laser  Division.  He  discussed  their  work  to  find 
efficient  diode  laser  pumped  solid  state  lasers  in  the  1  micro  band.  He  went 
into  some  detail  concerning  their  interest  in  activating  laser  crystals  for 
broadband  absorption  of  laser  diode  array  emission,  such  as  those  discussed 
in  the  article  by  A.  A.  Kaminskii.  (Phys.  Stat.  Sol.  87,  1 1  (1985). 

Dr.  Stroud  then  gave  a  general  introduction  to  the  work  going  on  in  this 
group  at  the  University  of  Rochester  involving  laser  dynamics  in 
homogeneously  broadened  laser  systems.  He  briefly  surveyed  the  recent 
progress  in  understanding  spiking  and  chaotic  instabilities  in  homogeneously 
broadened  lasers. 

Karl  Koch  then  gave  a  presentation  in  the  instabilities  that  he  has  been 
studying  in  dye  lasers,  both  intrinsic  instabilities  and  the  instabilities  that  are 
induced  by  external  modulation.  He  pointed  out  that  in  dye  lasers  and  other 
tunable  lasers  of  the  phonon  assisted  variety  are  not  well  described  by  simple 
ideal  4  level  models.  He  described  the  work  going  on  in  Rochester  and  in 
Professor  Haker's  group  in  Stutgard  on  improved  models  for  these  systems. 

Stephen  Chakmakjian  then  discussed  his  work  on  alexandrite,  another 
phonon  assisted  homogeneously  broadened  laser.  He  described  in  some 
detail  modulation  spectroscopic  techniques  developed  in  Rochester  that 
allow  one  to  easily  measure  the  various  lifetimes  involved  in  the  internal 
population  dynamics  in  these  lasers. 

Dr  Stroud  then  summarized  the  work  and  described  some  of  the  projects 
that  were  planned  or  underway  in  this  area  in  Rochester.  This  was  followed 
by  a  general  discussion  on  the  ways  that  the  modulation  techniques  might  be 
usefully  implemented  to  help  characterize  the  activated  crystals.  Dr.  Stroud 
said  that  he  was  not  really  set  up  in  hts  laboratory  at  present  to  systematically 
go  through  a  long  list  of  materials  and  characterize  them,  but  that  he  could 
easily  look  at  one  or  two.  Another  possibility  is  that  one  or  two  of  the 
students  from  Rochester  might  spend  a  few  days  at  NVEOC  to  help  in  setting 
up  an  experiment  there  to  make  such  measurements. 

Another  problem  in  the  Laser  Division  was  developing  a  simple  way  of  deter- 
mining  whether  a  potential  laser  material  is  homogeneously  or 
inhomogeneously  broadened.  Dr.  Sroud  discussed  some  of  the  usual 
techniques,  and  promised  to  send  Dr.  Pir.to  a  writeup  summarizing  the 
problem 


CENTER  FOR  OPTO-ELECTRONIC  SYSTEMS  RESEARCH 
INTRODUCTION  AND  OVERVIEW 


LASERS  ARE  NOISIER  THAN  THEY  SHOULD  BE 


200  ps/cm 


5  ps/cm 


1  ps/cm 


100  ns/cm 


Spiking  in  Ruby  Laser  output 


ARE  LASERS  NATURALLY  UNSTABLE? 

•  Described  by  NONLINEAR  equations. 

•  General  solutions  do  not  exist. 

•  Behavior  may  be  qualitatively  different 
for  varying  parameters  and  initial 
conditions. 

•  Progress  has  been  made  in  study  of 
driven  nonlinear  systems. 

•  System  may  be  stable  in  some  regions 
of  parameters  by  display 
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for  other  values  of  parameters. 


Homogeneously  Broadened  Lasers 


Textbook: 

"An  ideal  homogeneously  broadened  laser  will  operate 
in  one  mode.” 

Engineer: 

"Nonsense!” 

Question: 

Are  broadband  tunable  lasers  ideal  homogeneously 
broadened  systems? 
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HIGH-Q  RING  DYE  LASER 


5cm  radius 


< - 5cm 


Low  loss 

1  Watt  pump  produces  50  Watts  circulating  dye  laser 
power 

10  \xm  focal  spot  generates  50  MWatts  /cm2  in  the  gain 
region 

Pinhole  ensures  single  transverse  mode 


WAVELENGTH 


TWO-FREQUENCY  INSTABILITY 

CHARACTERISTICS 

Low  threshold  <<  9  times  above  threshold 

Two-frequency  operation  whose  splitting  increases  as 
the  square  root  of  the  intensity. 

Discontinuous  jumps  in  output  intensity  and  associated 
hysterisis  loops  as  a  function  of  pump  power. 
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EFFECTS  OF  DETUNING 
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UNIDIRECTIONAL  OPERATION 


5cm  rodius 


Bidirectional  Unidirectional 

The  two  frequency  instability  is  still  present  in  the 
unidirectional  laser. 

The  quartz  crystal  plate  acts  as  an  etalon  and  imposes  a 
coarse  mode  structure  on  the  spectrum. 
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EFFECTS  OF  FREQUENCY 
DEPENDENT  LOSS 


CAVITY  POWER  CAVITY  POWER 


340A  SPLITTING 
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TWO-FREQUENCY  ANALYSIS 


We  analyze  the  atomic  response  to  a  field  made  up  of  two 
frequencies  of  equal  amplitude.  The  two  components  are 
separated  by  28co.  8co  is  referred  to  as  the  modulation 
frequency. 


The  saturated  gain  experienced  by  the  two-frequency 
field  exhibits  a  dip  near  zero  modulation  frequency. 


Maximum  gain  occurs  for  a  nonzero  modulation 
frequency. 
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For  large  field  amplitudes  the  dip  in  the  saturated  gai 
develops  structure.  This  structure  shows  subharmoni 
resonances  of  the  Rabi  frequency. 
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Ix  =  2fc*7;ra<r*a  180  (intensity) 


BAND  MODEL  EQUATIONS  OF  MOTION 


BROADBAND  OUTPUT 


Power  Spectra 
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Power  Spectra 
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Driving  nonlinear  systems  with  two  incommensurate 
frequencies  can  produce  broadba:  d  output 


Inversion  In  |F(w)| 


RESULTS  OF  NUMERICAL  CALCULATION 
COMMENSURATE  MODULATION  FREQUENCIES. 
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RESULTS  OF  NUMERICAL  CALCULATION 
INCOMMENSURATE  MODULATION 
FREQUENCIES 
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Dye  Laser  Output  vs.  Argon  Laser  Input: 
modulation  above  threshold 


EXPERIMENTAL  SETUP 


logI0(argon  laser 
power  spectrum) 


COMMENSURATE  MODULATION 
FREQUENCIES 


Input  spectrum 
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Input  time  series 
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CENTER  FOR  OPTO-ELECTRONIC  SYSTEMS  RESEARCH 
MODULATION  TECHNIQUES:  ALEXANDRITE  LASERS 


Modulation  Techniques  in  Homogeneous 
Systems  to  Determine  Decay  Rates 


hole  burning 


inhomogeneous  broadening 


homogeneous  broadening 


pump-probe  techniques 

AM  modulation  techniques 

AM  sensitivity  (lock-in  detection) 
FM  noise  insensitive 
single  laser  experiment 


Applications  of  AM  Spectroscopy 


probe  gain,  or  absorption  in  presence  of 
saturating  field 


single  mode  stability  analysis 


measure  "hidden”  non-radiative  decay 
rates 


sensitivity  of  gain  transition  to  pump 
noise 


Three  Level  Systems 


Fluorescent  Decay  Rate 
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Modulation  Measurement 


Four  Level  Systems 
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fluorescent  emission  rate  measurement 
cannot  determine  if  non-radiative  decay 
from  level  1  is  slow 


modulation  measurement  determines  the 
population  cycling  rate 


slow  decay  from  level  1  may  cause 
instabilities 


Single  Laser  AM  Experiments 
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Results  of  Single-Laser  AM  Experiments 


a  hole  in  the  absorption  spectrum  of  the 
probe  fields  is  observed  centered  at  DC 


the  hole  has  width  1/Ti  for  weak  DC  field 
strength 


hole  broadens  as  1  +  I  /  ISat 
discrepancy  with  known  1/Ti  rate  in  ruby 

alexandrite  reveals  anti-hole  due  to 
excited  state  absorption 


Alexandrite  Level  Structure 


can  operate  as  a  three  level  laser  similar 
to  ruby  J ,  b  mi 

can  operate  as  a  four  level  laser  similar  to 
dyes,  a  vibronic  laser  u  (K  i 
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Rate  Equation  Approximation 

consider  a  two  level  system 

optical  Bloch  equations: 
u  =  -u/T2  -  Av 
v  =  Au  -  v/T2  +  xE(t)w 

W  =  -(w-WeqVTj  -  KE(t)v 
collisional  broadening,  1/T2  >  >  1/TL 
rate  equation  limit 


Reduced  Rate  Equations 
Three  Level  System 
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Four  Level  System 
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Multi-Level  Two  Laser  AM  Spectroscopy 


Experiments 

modulate  pump  to  study  pump  noise 

modulate  laser  tuned  to  inverted 
transition  to  study  single  mode  stability 


Magnitude  and  Angle 


Variation  of  Half-Widths  with  Intensity 
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Abstract 


A  series  of  resonances  has  been  observed  in  the  absorption  of  a  100% 
amplitude-modulated  laser  beam  by  an  optically  pumped  sodium  atomic  beam. 
Resonances  were  observed  when  the  modulation  frequency  was  equal  to  the 
first,  second,  and  third  subharmonic  of  the  Rabi  frequency.  The  experimental 
results  are  compared  with  theory  and  the  implications  of  these  parametric 
resonances  to  laser  instabilities  are  discussed. 
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INTRODUCTION 

Researchers  in  the  fields  of  modulation  spectroscopy,  optical  bistability  and  laser 
instability  have  become  interested  in  the  interaction  of  intense,  strongly  modulated  fields 
with  atomic  systems 1'7.  A  100%  amplitude-modulated  (AM)  field  is  a  special  case  of  strong 
modulation  since  the  energy  at  the  carrier  frequency  has  been  completely  transferred  into 
the  modulation  sidebands.  In  this  paper  we  study  the  interaction  of  a  100%  AM  laser  field 
with  an  ensemble  of  homogeneously  broadened  two-level  atoms. 

Under  the  influence  of  a  strong  resonant  field  the  atomic  variables  undergo  Rabi 
oscillations.  If,  in  addition  to  the  strong  pump  field,  a  weak  probe  field  is  applied,  the  probe 
field  will  see  a  complicated  structure  in  its  absorption  spectrum  due  to  the  presence  of  the 
strong  field.  When  it  is  detuned  from  the  strong  field  by  the  strong-field  Rabi  frequency  the 
probe  field  experiences  a  resonance  in  its  absorption  spectrum.8  10  These  results  have  been 
experimentally  confirmed 1113 .  Sargent  et  al. 14  have  shown  how  modulation  techniques 
can  be  employed  to  probe  the  saturated  absorption  of  these  strongly  driven  atoms. 
Modulation  techniques  alleviate  the  need  for  a  second  laser.  Instead,  the  strong  field  is 
weakly  modulated  at  a  tunable  frequency  and  the  modulation  sidebands  probe  the 
saturated  absorption.  Hillman  et  a/.15  and  Kramer  et  al.16  showed  that,  the  behavior  seen  in 
the  absorption  spectrum  of  the  probe  field  is  a  result  of  the  atomic  response  to  the 
modulation  tone  in  the  field  intensity.  In  their  analysis  they  treated  the  interaction  with 
the  strong  field  to  all  orders,  while  the  interaction  with  the  sidebands  was  treated 
perturbatively.  In  the  case  of  multiple  strong  fields  it  is  necessary  to  treat  all  the  field 
components  and  the  combination  tones  to  all  orders. 

Although  the  treatment  of  strong-modulation  interactions  is  by  no  means  complete, 
there  has  been  experimental  work  in  this  area.17,18  These  strong-modulation  experiments 
yield  different  results  from  those  seen  in  weak-modulation  experiments.  The  atomic 
variables  exhibit  resonant  behavior  when  the  modulation  frequency  is  approximately 
equal  to  the  Rabi  frequency  or  any  subharmonic  of  the  Rabi  frequency.  Bonch-Breuvich  et 
al.xl  performed  an  experiment  with  two  strong  radio-frequency  fields  tuned  to  a  Zeeman 
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resonance  in  cadmium.  One  radio-frequency  field  was  held  at  a  constant  frequency  as  the 
other  field  was  tuned  about  the  resonance  to  obtain  an  absorption  spectrum.  The 
absorption  spectrum  of  the  rf  fields  exhibited  several  subharmonic  resonances.  Thomann18 
performed  a  three-field  experiment  at  the  sodium  D2  line  resonance.  The  field  was 
produced  by  strongly  modulating  a  laser  beam;  the  Rabi  frequency  was  varied  while  the 
modulation  frequency  was  held  fixed.  Although  a  center  frequency,  the  carrier,  was 
present  in  this  experiment  the  physics  is  similar  to  the  two-field  experiments.  At  this  point 
it  is  helpful  to  review  the  theoretical  predictions  for  the  problem.4  6,19  22 
THEORETICAL  REVIEW 

The  response  of  an  atomic  system  can  be  characterized,  in  part,  by  the  rat-  at  which  it 
absorbs,  or  scatters,  energy  from  a  resonant  laser  beam.  We  employ  a  calculation 
developed  by  Hillman  et  al 6  to  calculate  the  rate  of  absorption  from  the  optical  Bloch 
equations 


u  =  -  u/T2  -  Au, 

(la) 

v  =  Au  -  vIT2  +  )aE(t)w, 

(lb) 

w  —  -(w  +  l)/Tj  -Kj£(t)u, 

(lc) 

where  A  is  the  detuning,  1/T2  and  i/T;  are  the  transverse  and  longitudinal  decay  rates 
respectively.  The  dipole  coupling  constant,  k,  is  given  by 

k  =  2  dJh  (2) 

where  d  is  the  dipole-moment  matrix  element  and  h  is  Planck’s  constant.  We  are  interested 
in  the  response  of  the  atom  to  a  100%  AM  excitation.  In  this  case  the  electric  field 
amplitude  can  be  written  as 

E(t)  =  2E]cos(8ut) ,  (3) 

where  5a  is  the  modulation  frequency  and  the  amplitude  of  the  electric  field  is  Er 

According  to  Floquet’s  theorem,  the  stationary-state  response  of  the  atomic  variables 
can  be  expanded  in  a  Fourier  series  of  the  modulation  frequency  5a.  The  atomic 
polarization  components  and  the  atomic  inversion  can  be  written 
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The  n^h  harmonic  componei.  cs  of  the  in-phase  and  in-quadrature  dipole  moments  and  of  the 
atomic  inversion  are  denoted  by  un,  vn,  and  wn,  respectively.  We  substitute  Eqs.  (4)  into  the 
Bloch  equations,  Eqs.  (1),  to  obtain  a  set  of  recurrence  relations  between  the  harmonic 
components  of  the  atomic  variables.  The  ratio  of  the  first  harmonic  component  of  the  in¬ 
quadrature  dipole  moment  to  the  zeroth  harmonic  (time-averaged)  component  of  the 
inversion  is  expressed  as  a  continued  fraction 


(5) 


B, + 
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where  the  single  field  intensity  is 


and  the  coefficients  denoted  by  B  are 
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for  odd  n,  and 
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B  =  1  4-  inScoT.  , 
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(6) 


(7a) 


(7b) 


for  even  n. 

We  can  solve  for  the  time-averaged  component  of  the  atomic  inversion,  w0,  in  terms  of 
the  continued  fraction,  and  then  write  the  excited  state  population 

=  2yMy%)  (8) 

p'  1  +2IxRe\uy/ 
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We  are  ultimately  interested  in  the  absorption  of  the  AM  field.  The  average  rate  at  which  a 
two-level  atom  scatters  light  in  the  form  of  resonance  fluorescence  is  proportional  to  the 
time-averaged  absorption.  The  fluorescent  emission  rate  is  proportional  to  the  excited 
state  population.  Consequently,  for  the  case  of  a  100%  AM  excitation  of  a  two-level  atom, 
the  rate  of  absorption  of  energy  is  proportional  to  the  time-averaged  component  of  the 
excited  state  population. 

We  plot  the  results  of  these  calculations  in  Figs.  1  and  2.  In  Fig.  1  we  plot  the  time- 
averaged  excited-state  population  for  a  fixed  Rabi  frequency  as  a  function  of  the 
modulation  frequency,  ScoTr  This  curve  is  quite  complicated  and  it  shows  several 
resonances  in  the  absorption  spectrum  of  the  fields.  To  understand  these  peaks  it  is  helpful 
to  solve  the  problem  in  the  absense  of  damping.  In  this  case  one  obtains  an  analytic 
solution  in  which  the  time-averaged  response  of  the  population  can  be  written  in  terms  of 
J0(2kE1/Sq).  When  we  numerically  pick  off  the  position  of  each  peak  from  Fig.  1,  we  find 
that  the  peaks  occur  whenever  the  factor  J0(2kEj/8co)  is  equal  to  zero.  The  outermost  peaks 
correspond  to  the  first  zero  of  JQ.  The  next  peaks  at  smaller  values  of  8qT2  correspond  to  the 
second  zero  of  J0  and  are  the  first  subharmonic  resonances.  In  the  limit  of  large  arguments, 
x,  the  zeros  of  J0(x)  become  equally  spaced  versus  x.  It  is  in  this  limit  that  the  resonances 
are  best  described  as  subharmonic  resonances.  In  Fig.  2  the  time-averaged  population  of 
the  excited  state  is  plotted  versus  the  time-averaged  dimensionless  intensity,  2(kE1)cT1T„, 
for  a  fixed  modulation  frequency  of  three  atomic  linewidths.  This  is  the  form  of  our  raw 
data.  The  resonances  appear  as  enhancements  in  the  absorption  embedded  in  the  usual 
single-field  saturation  curve.  As  the  intensity  is  increased  from  zero  we  approach  the  first- 
order  resonance.  The  second  resonance,  occuring  at  a  higher  intensity,  is  the  first 
subharmonic  resonance.  The  order  of  the  subharmonic  resonances  increases  as  we  look  to 
higher  dimensionless  intensities.  These  resonances  occur  only  when  the  modulation- 
interaction  is  strong,  meaning  the  Rabi  frequency,  kEp  exceeds  the  natural  linewidth  of 
the  transition,  and  the  modulation  depth  is  large,  as  in  the  case  of  100%  AM. 
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The  dynamic  behavior  of  the  atomic  variables  also  reflects  the  subharmonic  behavior 
of  the  100%  AM  interaction.  In  Figs.  3(a)-3(d),  we  show  the  phase  plots  formed  by  the  in¬ 
quadrature  atomic  polarization  and  the  inversion  at  each  of  the  first  four  resonances  (on- 
resonance  excitation  is  used  for  the  phase  plots  so  that  only  the  in-quadrature  polarization 
is  driven).  Each  figure  depicts  the  trajectory  followed  by  the  atomic  variables  for  a 
complete  period  of  the  modulation.  Fig.  3(a)  depicts  the  period-one  behavior  which  occurs 
at  the  first  resonance.  Fig.  3(b)  shows  that  a  second  cycle  in  the  trajectory  occurs  at  the 
second  resonance,  which  is  the  first  sub-harmonic  resonance.  Figs.  3(c)  and  3(d)  show  the 
behavior  at  the  third  and  fourth  resonances  respectively.  We  can  compare  this  system  to  an 
anharmonic  oscillator  being  driven  every  nth  cycle  of  its  natural  frequency.  The  oscillator 
can  be  driven  effectively  by  the  fundamental  frequency  or  any  subharmonic  of  the 
fundamental.  This  subharmonic  driving  is  effective  since  the  system  responds  at  higher- 
harmonics  of  the  driving  frequency  when  the  driving  force  is  large  enough  to  induce  a 
nonlinearity. 

EXPERIMENTAL  SETUP 

The  experimental  apparatus  we  used  to  study  the  absorption  spectrum  is  shown  in 
Fig.  4.  The  absorption  is  determined  by  measuring  the  total  fluorescence  from  a  small 
portion  of  the  interaction  region.  We  imaged  the  fluorescence  from  the  interaction  region 
onto  a  pinhole  placed  in  front  of  a  photo-multiplier  tube.  Spatially  filtering  the 
fluorescence  signal  limits  the  signal  to  that  from  a  small  region  near  the  center  of  the  laser 
beam  where  the  intensity  is  constant  to  within  10%.  This  allows  us  to  study  atoms  with 
equal  Rabi  frequencies.  The  detectors  were  dc-coupled  to  measure  the  time-averaged 
fluorescence. 

Sodium  was  chosen  as  the  atomic  medium  because  it  has  a  large  oscillator  strength 
and  is  readily  made  into  a  two-level  atomic  system.  To  obtain  a  two-level  atomic  system  we 
resonantly  excited  the  3s  to  3p  transition  with  circularly-polarized  light.  A  two-level 
system  results  from  the  pumping  of  population  into  the  aligned  magnetic  sublevels23.  The 
atomic  beam  was  collimated  to  provide  a  divergence  angle  of  1  milliradian  with  an 
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associated  Doppler  width  of  less  than  2  MHz  (this  is  20%  of  the  natural  line  width  of  the 
sodium  transition).  The  atoms  were  excited  by  a  circulariy-poiarized  beam  whose  angular 
divergence  was  less  than  1  milliradian  (this  angular  spread  accounts  for  no  more  than  1 
MHz  Doppler  broadening).  The  laser  field  was  the  output  of  a  frequency  stabilized 
rhodamine  6G  dye  laser,  whose  full-width-frequency  jitter  was  less  than  1  MHz.  The 
angular  divergence  of  the  laser  and  atomic  beams  together  with  the  frequency  jitter  of  the 
laser  gave  us  a  total  systematic  broadening  of  4  MHz,  or  40%  of  the  natural  linewidth.  This 
systematic  linewidth  was  verified  by  measuring  the  full  spectral  width  of  the  atomic 
transition  to  be  14  MHz  (the  natural  linewidth  is  10  MHz). 

We  produced  a  100%  AM  laser  field  with  a  modulation  frequency,  8co/2n,  of  15MHz, 
which  is  three  times  the  transverse  relaxation  rate  of  the  sodium  D2  line.  The  modulation 
was  produced  by  passing  the  laser  through  a  high  efficiency  acousto-optic  modulator 
(AOM).  A  portion  of  the  beam  was  up-shifted  by  the  30  MHz  drive-frequency  of  the  AOM. 
The  diffraction  efficiency  was  adjusted  so  that  the  power  diffracted  into  the  first-order 
beam  was  equal  to  that  left  in  the  zeroth-order  undiffracted  beam.  We,  then  recombined 
the  two  orders  to  create  the  bichromatic  field.  The  two  beams  were  aligned  with 
interferometric  precision  in  order  that  a  strong  beat  note  of  30  MHz  was  clearly  detected  by 
a  photo-diode  measuring  the  far-field  intensity.  After  recombination,  the  field  state  can  be 
described  in  terms  of  a  earner,  at  the  mean  frequency  of  the  two  frequency  components, 
that  is  100%  amplitude-modulated  at  one  half  the  drive  frequency  of  the  AOM.  Since  the 
recombination  geometry  is  different  for  each  modulation  frequency  it  was  inconvenient  for 
us  to  vary  the  modulation  frequency  in  search  of  the  resonances.  Instead,  we  held  the 
modulation  frequency  fixed  and  varied  the  Rabi  frequency  by  sweeping  the  intensity  of  the 
laser.  This  was  done  with  a  set  of  linear  polarizers  and  an  electro-optic  cell  driven  by  our 
computer. 

At  each  power  setting  we  recorded  the  incident  laser  intensity  along  with  the 
intensity  of  the  fluorescence  emitted  by  the  atoms.  Each  data  point  was  averaged  100 
times  over  several  milliseconds  to  integrate  out  any  fast  intensity  fluctuations  of  the  laser 
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field.  All  data  was  recorded  with  a  14  bit  analog  to  digital  converter  on  a  PDP-11/23 
microprocessor.  The  40  db  dynamic  range  provided  by  the  digitization  is  necessary  to 
match  the  dynamic  range  of  our  data. 

EXPERIMENTAL  RESULTS 

We  collected  data  for  the  AM  field  on-resonance  and  for  a  detuning  of  three 
linewidths.  The  modulation  frequency  was  equal  to  three  atomic  linewidths  for  all  the  data 
we  present  in  this  paper  (SqT2  =  3).  In  Fig.  5  we  plot  the  fluorescent  intensity  versus 
incident  time-averaged  intensity  for  on-resonance  excitation.  The  resonances  appear  as 
bumps  on  the  saturated  absorption  curve.  The  solid  curve  shows  a  theoretical  fit  to  our 
data.  A  two-level  theory  is  averaged  over  detuning  and  modulation  frequency  to  model  the 
very  slight  imperfect  recombination  of  the  two  beams.  When  the  two  laser  beams  exiting 
the  AOM  are  not  perfectly  aligned  (see  Fig.  4),  the  atoms  see  a  moving  interference  pattern. 
If  this  pattern  moves  along  the  atomic-beam  axis,  the  spread  in  atomic  velocities  causes  the 
atoms  to  experience  a  spread  in  modulation  frequencies. 

Since  we  are  mainly  interested  in  the  resonances,  we  can  divide  our  data  by  a  single¬ 
frequency  absorption  curve  to  obtain  a  normalized  atomic  response  to  the  two-field 
excitation.  This  normalization  flattens  out  the  curve  so  that  the  resonances  may  be  seen  on 
a  single  expanded  graph.  In  Fig.  6  we  plot  the  normalized  fluorescent  intensity  versus  the 
parameter  V^kLjTj.  The  parameter  V^TkEj  is  the  root-mean-square  (rms)  Rabi 
frequency.  A  distinct  enhancement  in  the  absorption  of  the  fields  occurs  when  the  rms  Rabi 
frequency  is  equal  to  n(28co),  where  n  is  the  number  of  the  resonance  and  28co  is  the 
frequency  difference  between  the  two  field  modes.  We  can  observe  the  first  order  resonance 
and  the  first  two  subharmonic  resonances  with  resonant  excitation.  With  the  fields 
detuned  by  three  natural  linewidths  we  have  observed  up  to  three  subharmonic  resonances. 
In  Fig.  7  we  have  plotted  the  fluorescent  intensity  versus  the  rms  Rabi  frequency  for 
resonant  and  detuned  excitation  (  for  the  detuned  case  AT2  =  3).  The  resonances  occur  at 
lower  intensities  when  the  modulated  field  is  detuned  since  the  generalized  Rabi  frequency 
is  larger  for  the  same  amount  of  incident  laser  power.  We  do  not  have  an  analytical 
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expression  for  the  position  of  the  resonance  peaks.  However,  we  find  that  the  empirical 
formula 

1  (9) 

(VFkjEj)2  +  A21*  =  n  (2Sa>), 

describes  the  postion  of  the  resonances  quite  well  for  both  detuned  and  resonant  excitation. 
We  can  interpret  the  quantity  on  the  left  as  the  generalized  Rabi  frequency.  A  resonance 
occurs  whenever  the  generalized  Rabi  frequency  is  equal  to  an  integer,  n,  times  the 
frequency  separation,  26co,  between  the  two  fields.  This  formula  for  the  resonances  agrees 
with  weak  modulation  experiments11’13. 

It  is  instructive  to  look  at  the  positions  of  the  resonances  for  the  detuned  and  resonant 
modulated  excitations.  We  used  the  computer  to  determine  the  position  of  the  resonance 
peaks  for  the  normalized  data.  The  positions  of  these  peaks  for  detuned  and  resonant 
excitation  are  plotted  on  the  same  graph  in  Fig.  8.  The  position  of  the  nth  resonance  for  the 
detuned  case  appears  at  a  lower  intensity  than  the  corresponding  resonance  for  the  case  of 
zero  detuning.  The  solid  lines  represent  the  theoretical  predictions  for  the  positions  of  the 
resonances  according  to  Eq.  (9). 

CONCLUSIONS 

We  have  presented  both  theory  and  experiment  for  the  absorption  of  a  100%  AM  field 
by  a  closed  two-level  atomic  system.  The  interaction  reveals  resonances  at  subharmonics  of 
the  Rabi  frequency.  We  observed  the  first  three  resonances  with  resonant  excitation,  and 
we  observed  the  first  four  resonances  with  detuned  excitation.  For  resonant  excitation,  the 
absorption  of  the  AM  field  is  enhanced  whenever  the  ratio  2kE1/8g)  is  equal  to  a  zero  of  the 
zeroth  order  Bessel  function.  In  this  paper  we  have  presented  an  empirical  relation  for  the 
occurence  of  the  resonances  for  both  resonant  and  detuned  excitation.  The  data  agrees  with 
this  empirical  prediction. 

Our  data  shows  the  absorption  of  the  bichromatic  field  is  enhanced  at  the  first-order 
resonance  and  at  each  of  the  subharmonic  resonances.  This  result  differs  from  that 
obtained  by  Thomann.18  The  results  of  that  experiment  showed  a  decrease  in  the 
absorption  at  each  resonance.  In  that  experiment  the  field  modulation  index  was  no  larger 
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than  0.65  so  there  was  always  a  significant  field  component  at  the  carrier  frequency.  When 
the  carrier  frequency  is  present  it  can  beat  with  each  of  the  modulation  sidebands  to 
produce  overtones,  at  harmonics  of  the  modulation  frequency,  in  the  atomic  response.  This 
beat  frequency  is  one  half  the  beat  frequency  arising  between  the  two  modulation 
sidebands.  The  harmonic  overtones  of  these  two  sets  of  beat  notes  can  destructively 
interfere  in  a  three  field  experiment.  This  interference  causes  the  absorption  to  diminish  at 
each  resonance. 

We  have  demonstrated  the  complicated  behavior  of  a  purely  bichromatic  field 
interacting  with  a  closed  two-level  atomic  system.  Subharmonic  resonances  occur  due  to 
the  coherent  nonlinear  interaction  of  a  multi-  frequency  field  with  an  isolated  atomic 
resonance.  It  is  exactly  this  type  of  nonlinear  interaction  that  determines  the  competition 
or  cooperation  between  cavity  modes  in  a  laser  gain  medium.  The  same  set  of  subharmonic 
resonances  occur  in  a  saturated  two-level  amplifier  interacting  with  a  modulated  laser 
field.  This  provides  a  mechanism  for  multimode  instabilities  in  a  homogeneously 
broadened  laser.  Instabilities  of  this  type  have  been  observed.5 
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Figure  Captions 


Fig.  1  The  time- averaged  excited  state  population  as  a  function  of  modulation 
frequency,  SC0T2.  The  Rabi  frequency  is  held  fixed,  kEjT2  =  10,  ^hile  the 
modulation  frequency,  ScjT2,  is  varied  to  observe  the  subharmonic  resonances. 
The  detuning  is  zero,  AT2  =  0. 

Fig.  2  The  time-averaged  excited-state  population  as  a  function  of  intensity.  The 
modulation  frequency  is  held  fixed,  8coT2  =  3,  while  the  time-averaged 
dimensionless  intensity,  2(kE1)2T1T2,  is  varied  to  observe  the  resonances.  The 
detuning  is  zero,  AT3  =  0. 

Fig.  3  The  trajectories  of  the  atomic  inversion  and  the  in-quadrature  polarization.  We 

plot  th®  atomic  inversion  versus  the  atomic  polarization  at  the  first  four 
resonances.  The  modulation  frequency  for  each  of  the  figures  is  8coT2  =  3.  In  Fig. 
3(a)  the  Rabi  frequency  corresponds  to  the  first  resonance.  The  next  three 
figures  show  the  behavior  at  the  first  three  subharmonic  resonances  (b)  n  =  2,  (c) 
n  =  3,  and  (d)  n  =  4.  The  detuning  is  zero,  AT3  =  0. 

Fig.  4  The  experimental  apparatus.  The  diagram  is  labeled  with  the  following 
notation:  SDL  -  stabilized  dye  laser;  AOM  -  acousto-optic  modulator;  M  - 
mirrors;  BS  -  beam  splitter;  LP  -  linear  polarizer;  EOM  -  eletro-optic  modulator; 
A/4  -  quarter-wave  plate;  Na  -  sodium  oven;  AB  -  atomic  beam;  PH  -  pin-hole; 
PMT  -  photo-multiplier  tube;  PD  -  photo-diode;  PDP-11/23-  micro-computer. 
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Fig.  5 


Fig.  6 


Fig.  7 


Fig.  8 


The  time-averaged  fluorescent  intensity  versus  time-averaged  incident 
intensity.  The  raw  data  is  plotted  with  dots  and  the  solid  curve  represents  a  best 
theoretical  fit.  The  modulation  frequency  is  fixed,  8uT2  =  3,  and  the  detuning  is 
zero,  AT2  =  0. 

: 

The  normalized  time-averaged  fluorescent  intensity  versus  rms  Rabi  frequency 
for  resonant  excitation.  The  fluorescence  data  from  Fig.  5  was  normalized  by  a 
single  field  saturation  curve  to  flatten  out  the  resonant  peaks.  The  Rabi 
resonance  and  the  first  two  subharmonic  resonances  are  shown. 

Normalized  time-averaged  fluorescent  intensity  versus  rms  Rabi  frequency  for 
resonant  and  detuned  excitation.  The  dotted  curve  represents  the  data  for  the 
resonant  case,  AT2  =  0.  The  dashed  line  represents  data  for  the  detuned  case, 
AT2  =  3.  The  resonant  peaks  occur  at  lower  Rabi  frequencies  for  the  detuned 
case.  The  modulation  frequency  is  8uT2  =  3. 

The  position  of  the  subharmonic  resonances  versus  resonance  number  for 
resonant  and  detuned  excitation.  The  rms  Rabi  frequency,  VJ'kEjTj, 
corresponding  to  the  n^  resonance  is  plotted  versus  n.  The  circles  represent  the 
data  for  resonant  excitation,  AT2  =  0.  The  squares  represent  the  data  for 
detuned  excitation,  AT2  =  3.  The  solid  line  represents  the  empirical  expression 
in  Eq.  (9)  for  resonant  excitation  AT2  =  0.  The  dashed  line  represents  Eq.  (9)  for 
AT2=  8to.  The  modulation  frequency  is  8uT2  =  3. 
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Abstract 

We  consider  the  interaction  of  two  100%  amplitude-modulated  fields 
with  a  two-level  system.  The  analysis  includes  a  discussion  of  the 
stability  of  an  amplitude-modulated  laser  field  to  the  growth  of 
subharmonic  amplitude-modulated  fields. 
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I.  Introduction 

Recently,  a  great  deal  of  interest  has  focused  on  instabilities  in  laser  systems1'2. 
Motivation  of  this  interest  is  provided  by  a  desire  to  understand  the  complicated 
dynamics  of  nonlinear  systems.  Optical  systems  provide  a  convenient^round  for 
rigorously  examining  these  dynamics.  The  laser  is  a  specific  example  of  a  rigorously 
describable  optical  system  containing  the  nonlinear  interaction  of  light  with  matter. 
Subharmonic  instabilities  have  been  observed  in  a  number  of  different  lasers3  5. 
Subharmonic  instabilities  have  become  especially  interesting,  since  a  sequence  of 
period  doublings  has  been  identified  as  a  possible  route  to  chaos. 

Instabilities  in  lasers  have  existed  since  the  first  demonstration  of  a  laser  by 
Maiman6.  Lasers  have  operated  in  a  number  of  states  since  that  first  demonstration. 
The  complications  introduced  by  pump  fluctuations,  spatial  hole  boning, 
inhomogeneous  broadening,  etc.,  have  enhanced  this  variety  of  available  states.  The 
continuous-wave  ring  dye  laser  is  a  homogeneously  broadened  laser  that  contains 
few  of  these  additional  complications,  yet  still  displays  a  rich  variety  of  operational 
states.  Hillman  et  al .5  observed  a  continuous-wave  dye  laser,  initially  operating  at  a 
single  frequency,  go  unstable  and  begin  oscillating  at  two  frequencies.  As  the  pump 
power  was  increased,  the  two  components  of  the  laser  spectrum  separated  with  a 
linear  dependence  on  the  electric  field  strength  of  the  laser.  When  the  separation  of 
the  spectral  components  reached  — 80A,  the  laser  encountered  another  instability 
threshold.  The  two-frequency  state  then  collapsed  to  another  two-frequency  state  of 
smaller  frequency  separation.  On  increasing  the  pump  power  still  further,  this  new 
two-frequency  state  behaved  similarly  to  the  initial  two- frequency  state. 

The  physical  mechanism,  for  these  new  states  of  dye  laser  operation,  is  not  well 
understood.  There  have  been  a  number  of  theoretical  treatments  of  a  multimode 
laser711.  Unfortunately,  none  of  them  has  provided  an  explanation  for  all  the 
features  of  the  experimental  results.  The  most  successful  theoretical  model3'1 1  takes 
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the  molecular  structure  of  the  laser-dye  molecule  into  account.  However,  the  actual 
physical  mechanism  responsible  for  the  instability  remains  to  be  seen. 

In  an  attempt  to  develop  physical  intuition  for  the  problem  of  a  nonlinear 
system  being  driven  by  a  two-component  field,  we  previously  examined^he  gain  for 
such  a  field  driving  a  two-level  atom8.  These  calculations  showed  that  the  gain  of  the 
two-component  field  displays  resonances  when  the  frequency  separation  of  the  two 
fields  is  equal  to  the  Rabi  frequency  or  a  subharmonic  of  the  Rabi  frequency. 

In  this  paper  we  consider  the  stability  of  a  bichromatic  laser  field.  Since  a 
general  self-consistent  solution  to  the  multimode  laser  equations  is  not  available, 
the  stability  analysis  must  be  carried  out  in  some  other  manner.  The  equations 
could  be  numerically  integrated;  however,  the  parameter  space  is  too  large  for  such  a 
solution  to  provide  intuition.  An  alternative  is  to  calculate  the  atomic  response  to 
the  laser  field  and  examine  the  gain  experienced  by  a  probe  field.  When  the  probe- 
field  gain  exceeds  the  laser-field  gain,  an  instability  will  occur.  In  this  paper  we  will 
discuss  the  stability  of  a  bichromatic  field  to  the  growth  of  a  bichromatic  probe  field 
using  just  such  an  alternative. 

The  analytic  description  of  the  interaction  of  four  field  components  wi  h  a 
resonant  medium  is  a  very  cumbersome  task.  The  nonlinear  response  of  the  medium 
produces  combination  tones  of  the  field  components  in  the  polarization  of  the 
medium.  When  the  frequencies  of  the  combination  tones  are  incommensurate  with 
one  another,  the  medium  can  exhibit  a  broadband  spectral  response52.  The  analytic 
description  of  such  a  problem  becomes  extremely  difficult.  However,  if  the 
frequencies  of  the  combination  tones  are  commensurate,  one  can  employ  Floquet’s 
theorem13  to  analytically  study  the  problem.  Neglecting  dispersion  in  the  cavity,  the 
mode  frequencies  are  all  commensurate  with  one  another.  The  mode  frequencies  are 
integer  multiples  of  the  inverse  of  the  cavity  round-trip  time.  Therefore,  the 
simplest  physical  situation  is  one  in  which  the  field  frequencies  are  commensurate 
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with  each  other.  The  most  elementary  four-component  field  is  one  in  which  each 
field  component  creates  a  degenerate  beat  frequency  with  the  nearest  other  field 
component.  We  will  examine  such  a  field  and  further  specify  the  problem  by 
considering  only  subharmonic  probe  fields.  These  probe  fields  are-especially 
interesting  as  they  give  us  insight  into  possible  physical  mechanisms  for 
subharmonic  instabilities.  This  same  analysis  could  be  carried  out  for  other 
subharmonic  fields,  such  as  the  -J-subharmonic  observed  in  period  doubling 
sequences.  The  techniques  used  in  the  analysis  of  such  a  field  would  be  identical, 
however  the  dimension  of  the  matrices  used  in  the  solution  would  be  larger.  The 
stability  of  bichromatic  fields  to  harmonic  bichromatic  probe  fields  will  be  the 
subject  of  a  future  paper14.  For  simplicity,  the  fields  are  symmetrically  placed  about 
the  atomic  resonance  frequency;  the  extension  of  the  problem  to  include  asymmetric 
field  states  is  straight  forward.  Figure  1  shows  the  spectrum  of  this  symmetric  field. 
2.  Outline  of  the  calculation 

The  general  expression  for  a  four-component  field  is  better  understood  by  first 
analyzing  the  single  100%  amplitude-modulated  (AM)  field.  The  expression  for  such 
a  field  is  written  as  a  two  component  field  or  a  single  component  with  a  sinusoidal 
amplitude. 

Ej(t)  =  2  E, (  cos(co  +  8cj)t  +  cos(co  -  5co)t )  (la) 

=  4  £  cosScot  coscat .  (lb) 

The  amplitude  of  each  frequency  component  is  equal  to  2£,. 

The  use  of  complex  phasor  diagrams  to  describe  the  phase  relationships 
between  many-component  fields  can  be  very  helpful.  The  complex  phasor  for  a  100% 
AM  field  is  a  vector  that  oscillates  along  a  line  through  the  origin.  The  amplitude  of 
the  field  determines  the  length  of  the  vector;  the  modulation  frequency  determines 
the  vector’s  rate  of  oscillation.  In  defining  the  field  above,  we  chose  an  arbitrary 
phase  factor  so  the  complex  phasor  lies  along  the  real  axis  of  the  complex  plane. 
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There  axe  two  independent  subharmonic  probe  fields  at  each  modulation 
frequency.  When  the  field  is  symmetric  about  the  atomic  resonance,  as  in  our  case, 
the  independent  probe  fields  are  in-phase  or  in-quadrature  to  the  strong  field;  these 
phasors  will  be  along  the  real  or  imaginary  axes  respectively.  The  expression  for  the 
100%  AM  perturbative  fields  will  be  similar  to  the  expression  above;  however,  in 
general  there  will  be  a  phase  difference  between  the  oscillations  of  the  perturbative 
fields  and  the  strong  field.  With  these  considerations  in  mind,  we  choose  the 
perturbation  to  the  100%  AM  field  for  our  calculations  to  be  of  the  form 

8E(t)  =  4  8£j'  cos^Scot  +.0')  coscot  -  4  S^"sin(-^8cot  +  0")  sincot .  (2) 

In  Eq.  (2),  the  amplitude  of  the  in-phase  (in- quadrature)  probe  field  is  Sf,'  (S^").  The 
phase  between  the  oscillations  of  the  strong  field  and  the  in-phase  (in-quadrature) 
subharmonic  probe  field  is  0'  (0").  The  total  field  is  the  sum  of  Eq.  (1)  and  Eq.  (2). 

The  optical  Bloch  equations  model  the  response  of  a  two-level  atom  coherently 
driven  by  an  optical  field15.  In  a  reference  frame  rotating  at  the  frequency,  co,  the 


optical  Bloch  equations  are 

\i  =  -Q"(t)  w  -  u/T2 ,  (3a) 

y  =  fi'(t)  w  -  V/T2 ,  (3b) 

=  Q"(t)  u  -  fi'(t)  v  -  (w  -  w^VTj  ,  (3c) 

in  which 

fl'(t)  =  2k.[8£,'  cos(-3-8ut  4-  0')  +  ^cosScot]  ,  (4a) 

fT(t)  =  2k  8E,"  sin(-j8cot  +  0")  ,  (4b) 


and  k  =  2d Jk.  The  dipole-moment  matrix  element  between  the  two  levels  is  d  and  h 
is  Planck’s  constant.  The  slowly  varying  amplitudes  of  the  in-phase  (in-quadrature) 
part  of  the  dipole  moment  is  u  (v).  The  inversion  of  the  two-level  system  is  w;  w  is 
the  value  the  inversion  decays  to  in  the  absence  of  a  coherent  field.  The  relaxation 
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times  of  the  inversion  and  polarization  are  T:  and  T2,  respectively.  The  field 
components  that  are  in-phase  (in-quadrature)  to  the  rotating  frame  are  £V(t)  (Q"(t)). 

Floquet’s  theorem  tells  us,  in  the  stationary  state  the  atomic  variables  will 
respond  at  harmonics  of  the  driving  field  modulation  frequency.  -So,  in  the 
stationary  state  we  express  the  atomic  variables  as  Fourier  series  in  the  modulation 
frequency. 


u  =  ^  exp($in&(j}t)  ,  v=  ^  ezp(\Ln5u>t)  ,  w  =  wnezp(}ui&ut)  (5) 


This  assumption  leads  to  the  following  recurrence  relations 

(1  +  £inScoT2)  un  =  (8CT  wn.1  +  8 Q"*  wn  +  1)T2  ,  (6a) 

(1  +  ^inScoT2)  vn  =  (5Q'  wa  l  -I-  SQ '*  wn+1  +  Q'  wn3  +  Q'  wn+3)T2 ,  (6b) 

(1  +  -JinScoT^  wo  =  un4  +  5i'2  r*  un  +  1  -f  6fi'vnl  -c  SQ'*  vn+1 

+  n'  V-3  +n'  Vn+3>T1  *  (6c) 

in  which 

n'  =  ,  (7a) 

6Q'  =  k8E,'  exp(iQ')  ,  (7b) 

8Q”  =  iKS^”  exp(i9")  ,  (7c) 

and  the  asterisk  denotes  a  complex  conjugate.  The  Kronecker  delta  function  is  5q  q. 

We  solve  for  u  and  v  in  the  first  two  recurrence  relations  and  eliminate  them 
from  the  third  relation.  The  third  relation  then  becomes  a  function  of  the  Fourier 
components  of  w  alone. 


aw  +  b  w  .  +  c  w  +  d„  w  +  e  w  +  f  w„  .  .  +  g  w  . ,  =  w  8  . 
in  which 

an  =  Ln.3  fi'2TjT2  , 

\  =  +  Ln3)Q'8Q'TJT2  , 

cD  =  (LnVSQ"2  +  8Q'2)  +  (Ln  +  1  +  Ld3)  £2' BQ-)TXT2  . 


(8) 

(9a) 

(9b) 

(9c) 


J 
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dn  =  (1  +  -JinScdTj  +  (Lfl+l  +  Ln  l)  ((SQ'|‘  +  (6fl'|2)  T,T2 

+  (Ln-3  +  fl'2T,T2)  , 

e„  =  (Ldt1  <S!T *2  +  sn--2)  +  (L„.,  +  L„.,) fi'  60”)  T,T2  , 

f„=(Vi  +  LQ*3)n'8n"T.T2  ■ 

g„  =  L.o  n'2T,T2 , 

and 

Ln  =  (1  +  -JdnScoTj)1  . 

This  recurrence  relation,  Eq.  (8),  can  be  rewritten  as  a  three-term  recurrence 

relation;  the  terms  are  products  of  matrixes  and  vectors. 

PDWn-2  +  QnWn  +  RoWn.2  =  Dn 
in  which 


(9  d) 
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We  solve  this  three-term  recurrence  relation  using  matrix  continued  fractions16. 
The  solution  is  given  in  Appendix  A.  Now  that  we  have  determined  the  atomic 
response,  we  calculate  the  gain  experienced  by  the  probe  fields. 

The  time-averaged  product  of  the  probe  field  and  the  time  derivative  of  the 
polarization  gives  the  gain  experienced  by  the  probe  field.  The  gain  of  the  in-phase 
subharmonic  probe  field  is  proportional  to  the  real  part  of  the  first  harmonic  of  v. 

d  S£,'  /  dz  =  a  Re{Vj  exp(-i0')}  /  2kT2  ( 12a) 

The  gain  of  the  in-quadrature  subharmonic  probe  field  is  proportional  to  the 
imaginary  part  of  the  first  harmonic  of  u. 

d  8E,"  /dz  =  a  Im^expC-iQ')}  /  2kT2  (12b) 

3.  Discussion  of  results 

In  Fig.  2  we  plot  the  in-phase  subharmonic  probe-field  gain  as  a  function  of 
modulation  frequency  for  various  strong-field  intensities.  At  low  intensities  the 
probe-field  gain  is  a  Lorentzian  in  modulation  frequency  with  a  width  of  1/T2.  As  the 
intensity  of  the  strong  field  is  increased,  the  gain  curve  develops  a  dip  near  zero 
modulation  frequency.  The  dip  has  a  width  of  1/Tr  The  presence  of  this  dip  implies 
that  the  probe  field  will  experience  the  greatest  gain  for  a  nonzero  modulation 
frequency.  In  Fig.  3  we  plot  the  in-quadrature  subharmonic  probe-field  gain  as  a 
function  of  modulation  frequency  for  various  strong-field  intensities.  Again  we  see 
that  at  low  intensities  the  probe-field  gain  is  a  Lorentzian  in  modulation  frequency 
with  a  width  of  1/T2.  As  the  strong-field  intensity  is  increased,  the  gain  curve 
develops  an  anti-dip  near  zero  modulation  frequency. 

The  calculations  have  shown  both  quadratures  of  the  probe-field  gain  to  be 
independent  of  9'  and  9".  In  Appendix  B  wc  show  analytically  that  the  in¬ 
quadrature  suhharmonic  probe-field  gain  is  independent  of  9".  The  independence  of 
subharmonic  gain  from  the  phase  difference  between  the  oscillations  of  the  probe 
field  and  the  oscillations  of  the  strong  field,  implies  that  the  subharmonic  component 
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of  the  polarization  is  directly  proportional  to  the  subharmonic  field.  In  contrast,  it 
can  be  shown  that  the  harmonic  probe-field  gain  exhibits  a  strong  dependence  on  the 
phase  difference  between  the  strong-field  oscillations  and  the  harmonic  probe-field 
oscillations14.  Also,  if  the  amplitude  of  the  probe  fields  is  increased,  a»  that  their 
interaction  can  no  longer  be  considered  only  to  first  order,  the  gain  of  the 
subharmonic  field  becomes  dependent  on  0.  The  independence  of  the  gain  from  the 
phase  implies  the  instabilities  predicted  may  be  symmetry-breaking  instabilities. 
The  phase  of  the  light  field  for  a  laser  at  the  first  laser  threshold  displays  similar 
behavior17. 

In  Figs.  4  and  5  we  show,  respectively,  the  in-phase  and  in-quadrature 
subharmonic  probe-field  gains  for  larger  strong-field  intensities.  In  both  cases, 
additional  structure  begins  to  develop  about  the  initial  feature.  The  additional 
maxima  and  minima  seem  to  move  out  with  a  parabolic  dependence  on  intensity. 

In  Figs.  6  and  7  we  plot  the  peak  positions  of  the  subharmonic  gain  curves  as  a 
function  of  modulation  frequency  and  amplitude  of  the  strong  field.  Here  we  can 
clearly  see  the  peaks  move  out  to  higher  frequencies  linearly  as  the  electric  field 
increases.  This  linear  dependence  suggests  a  connection  between  these  resonances 
and  the  Rabi  frequency.  In  contrast  to  the  two-field  problem,  where  the  position  of 
the  resonances  are  found  to  be  independent  of  the  damping  rates8,  the  position  of  the 
subharmonic  probe-field  resonances  do  display  a  dependence  on  damping.  The 
maxima  occuring  at  the  largest  modulation  frequency  varies  most  drastically  with 
variation  of  the  ratio  between  Tj  and  T2.  The  subsequent  peaks  vary  less  and  less 
sensitively.  The  physical  mechanism  responsible  for  this  behavior  is  not  well 
understood.  Some  of  the  features,  however,  are  similar  in  behavior  to  that  seen  in 
multiphoton  absorption18.  Subharmonic  resonances  such  as  these  are  also  present  in 
the  two-frequency  problem  and  have  also  been  predicted  in  other  multichromatic 
analyses18. 
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In  Fig.  8  we  show  the  strong-field  gain  as  a  function  of  modulation  frequency 
for  a  fixed  strong-field  amplitude.  In  reference  8  the  maxima  in  the  strong-field  gain 
curve  were  found  to  occur  where  the  zeroth  order  Bessel  function,  J0(2k5  /  8co)  is 
equal  to  zero.  For  large  Rabi  frequencies,  k£,  and  small  modulation  frequencies,  8oo, 
the  zeroes  of  the  Bessel  function  become  equally  spaced.  It  is  in  these  regimes  that 
the  resonances  become  subharmonics.  In  Fig.  9  we  plot  the  in-phase,  in-quadrature, 
and  strong-field  gains  as  a  function  of  modulation  frequency  for  the  same  parameters 
as  Fig.  8.  At  this  intensity  the  in-quadrature  probe-field  gain  is  always  positive;  it 
exceeds  the  strong-field  gain  for  all  modulation  frequencies,  except  in  the  vicinity  of 
the  absolute  maxima  of  the  strong-field  gain  curve.  The  in-phase  probe-field  gain 
takes  on  both  positive  and  negative  values  as  a  function  of  modulation  frequency.  It 
can  also  be  noted  that  the  oscillations  in  modulation  frequency  of  the  two 
subharmonic  gain  curves  are  180°  out  of  phase  with  the  oscillations  in  modulation 
frequency  of  the  strong-field  gain.  It  is  important  to  remember,  however,  the 
frequency  separation  of  the  probe  fields  is  one  third  the  modulation  frequency,  8co. 
This  means  the  presence  of  a  frequency-dependent  loss  influences  the  two  curves 
differently.  For  instance,  a  frequency-dependent  loss  that  is  larger  for  larger 
frequencies,  reduces  the  strong-field  curve  more  than  the  probe-field  curve.  This 
reduction  makes  the  probe-field  gain  larger  than  the  strong-field  gain.  If  the  strong 
field  is  the  field  in  a  laser,  this  reduction  allows  the  probe  field  sufficient  gain  to 
overcome  losses  and  begin  oscillation.  The  severity  of  the  frequency-dependent  loss 
determines  the  point  at  which  the  instability  will  occur.  The  experiments  of  Hillman 
et  al  b  and  Stroud  et  al ,19  confirm  these  conclusions. 

In  Fig.  10  and  Fig.  11  we  plot  a  series  of  in-phase  and  in-quadrature 
subharmonic  probe-field  gain  curves  for  fixed  strong-field  amplitude  as  a  function  of 
modulation  frequency  and  three  different  ratios  of  damping  times.  The  ratios  chosen 
are  indicative  of  collisional  dephasing  (T2  <*  Tj),  strong  collisions  (T2  =  Tj),  and 


radiative  broadening  (T2  =  2  Tj).  These  plots  indicate  the  shape  of  the  in-phase  gain 
curve  is  not  a  strong  function  of  the  ratio  of  damping  times.  However,  as  noted 
before,  the  positions  of  the  maxima  do  show  a  dependence  on  the  ratio  between 
and  T2.  The  in-quadrature  subharmonic  probe-field  gain  does  indicate  ajiependence 
on  the  ratio  of  Tx  and  T2  for  modulation  frequencies  less  than  1/Tr  However,  for 
modulation  frequencies  greater  than  the  linewidth,  1VT2,  the  curves  are  quite  similar 
in  form  and  only  show  a  weak  dependence  on  the  ratio  of  Tj  to  T2. 

Conclusions 

We  have  calculated  the  atomic  response  to  a  pair  of  100%  AM  fields.  We  used 
this  solution  to  calculate  the  gain  seen  by  a  100%  AM  subharmonic  probe  field.  The 
subharmonic  probe  field  sees  a  Lorenthan  lineshape  in  modulation  frequency  for  the 
limit  of  small  amplitude.  As  the  strong-field  amplitude  is  increased  it  modifies  the 
lineshape  experienced  by  the  subharmonic  probe  field.  The  lineshape  contains  a 
primary  resonance  and  subharmonic  resonances  of  this  primary  resonance.  The 
position  of  the  resonances  varies  linearly  with  the  Rabi  frequency  of  the  strong  field. 
The  location  of  the  resonances  also  depends  on  the  ratio  of  the  inversion  decay  rate  to 
the  dipole  decay  rate. 

We  have  also  observed  the  gain  of  the  subharmonic  probe  field  to  be 
independent  of  the  phase  between  the  oscillations  of  the  strong  field  and  the 
oscillations  of  the  subharmonic  probe  field.  This  independence  of  phase  implies  the 
subharmonic  response  of  the  atom  is  directly  proportional  to  the  subharmonic  field. 
It  also  predicts  the  instabilities  will  be  independent  of  this  phase  difference. 

The  calculations  also  reveal  that  the  100%  AM  strong  field  is  stable  to  the 
growth  of  subharmonic  fields  if  there  are  no  frequency-dependent  loss  mechanisms. 
However,  if  frequency-dependent  loss  mechanisms  are  present,  we  predict  the  100% 
AM  field  will  go  unstable.  The  point  at  which  the  field  goes  unstable  depends  on  the 


severity  of  the  frequency-dependent  loss.  This  conclusion  agrees  well  with  observed 
instabilities  in  cw  dye  lasers619. 
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Appendix  A 

In  this  appendix  we  outline  the  method  for  calculating  the  gain  of  the  subharmonic 
probe  fields. 

The  larger  harmonic  components  of  w  must  become  less  signifieant  if  the 
Fourier  series  in  Eq.  (5)  are  to  converge  and  the  rotating-wave  approximation  is  to 
remain  valid.  Therefore,  we  assume  the  harmonic  components  of  the  inversion,  wd, 
are  negligible  for  n>N.  With  this  assumption,  the  recurrence  relation,  Eq.  (10),  for 
n  =  N  can  be  written  as 


PnWn;  +  Q.,.W,  =  0  ,  (Al) 

which  implies 

WN  =  *Qn  lpNwN-2  •  (A2) 

The  recurrence  relation,  Eq.  (10),  for  n  =  N-2  can  now  be  written  as 

^N-2  ^N-4  +  (Qn-2  *  ^N-2  ^N-2  =  ®  *  ^3) 

which  implies 

WN2  =  -(Q^-Rn^Qn^Pn^P^W^  .  (A4) 

We  continue  in  such  a  fashion,  defining  the  larger  index  Fourier  amplitudes  in  terms 
of  the  lower  index  Fou  er  amplitudes,  until  we  reach  the  relation  for  n  =  2.  This 
relation  has  the  following  form. 


W,  =  -<Q2  -  R2(Q4  -  R»(...(QK  2  -  RN2Q..,‘PS)  ‘  PN  2) 1  PN.,...)-'  Ps)-'  P.) 1  P2  W0,  (A5) 

or  more  simply 


w2  =  mw0  , 


(A6) 


where  M  is  the  matrix  continued  fraction,  with  elements  m^.  This  expression,  Eq. 
(A6),  relates  the  harmonic  components  w4,  w6,  and  wg  to  the  harmonic  components 
wQ,  w2,  and  w  2.  Strictly  speaking,  the  matrix,  M,  is  a  function  of  N.  However,  N  is 
chosen  sufficiently  large  so  that  M  becomes  independent  of  N  to  the  desired  accuracy 
of  the  calculation.  This  is  our  criterion  of  convergence.  Since  the  inversion,  w,  is  a 


real  quantity,  the  Fourier  components  obey  the  relation 


We  now  consider  Eq.  (8)  for  the  cases  n  =  -4,  -2,  0, 2,  and  4. 
(A7)  we  write  these  five  relations  in  matrix  form  as 
SW  =  W0 

Where  the  vectors  W  and  W0  are 


W  = 


i_w4 


and 


0 

0 

w 

0 

0 


And  the  matrix  S  has  elements  s^  where 

*11  =  d-4  +  a-4m33*  +  ^-4m23*  +  C-4  mi3*  * 
*12  ~  e-4  a-4m32*  ^-4m22*  C-4  mi2*  ’ 
S13  =  f-4  +  a-4m31*  +  ^-4m21*  +  C-4mil*  ’ 
SH  =  £-4  » 

*15  =  0' 

*21  =  C-2  "^a-2m23*  +  ^-2mi3*’ 

S22  =  d2  +  a.2m22*  +  b2ml2*  , 

*23  =  e-2  a-2m2]*  ^-2mil*  ’ 

*24  =  ^-2  * 

*25  =  £-2  * 

*31  ~  +  a0  mi3*  * 

*32  —  C0  a0mi2*  * 


(A7) 

Using  Eqs.  (A6)  and 

(A8) 

(A9) 


*33  —  ^0  ^  aomil*  » 


(A10) 
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534  =  S32*  » 

535  =  S31*  » 

541  =  *2  » 

542  =  ^2  ’  ” 

543  =  c2  +  f2mu  +  g2m,1, 

544  =  ^  +  f2  ml2  +  g2  mjj  , 

545  =  e2  +  *2  mi3  +  &2  “23  • 

s61  =  0  , 

552  =  a4  ’ 

553  =  b4  +  e4mil  +  +  *4  “31  » 

854  =  C4  +  e4mi2  +  f4“a2  +  84m32« 

S55  =  d4  +  e4mi3  +  *4^23  +  * 

Equation  (A8)  represents  a  system  of  five  linear  equations  with  five  unknowns.  The 
solution  gives  us  expressions  for  w0,  w2,  w4,  and  their  complex  conjugates.  We  can 
use  these  expressions  to  find  the  other  Fourier  components  of  w.  From  Eq.  (6a)  and 
Eq.  (6b)  we  can  then  determine  the  Fourier  components  of  u  and  v  that  determine  the 
gain. 

Appendix  B 

In  this  appendix  we  show  that  the  gain  of  the  in-quadrature  subharmonic  probe  field 
is  independent  of  the  phase  difference,  6*,  between  the  oscillations  of  the 
subharmonic  probe  field  and  the  strong  100%  AM  field. 

The  solution  for  the  atomic  response  to  a  single  resonant  100%  AM  field  has 
been  found  in  reference  8.  If  we  linearize  the  Bloch  equations  about  this  solution,  the 
equation  of  motion  for  the  perturbation  to  the  slowly  varying  amplitude  of  the  in- 
phase  dipole  moment  becomes 
5u  =  -8(1  "(t)  wAM  -  8u/T2  . 


(Bl) 
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Where  Su  is  the  perturbation  to  the  slowly  varying  amplitude  of  the  in-phase  dipole 
moment.  The  perturbative  field  in-quadrature  to  the  100%  AM  field  is  8£!"(t).  The 
stationary  solution  for  the  inversion  to  the  Bloch  equations  driven  by  a  resonant 
100%  AM  field  is  wAM.  Recall  wAM  can  be  written  as  ~ 

(B2) 

WAM  ~  X,  Wn  ■ 

n  m  -«» 

If  8Q"(t)  =  2Sfl"  sin(^8ojt  +  0")  we  can  employ  Floquet’s  theorem  and  write  the 

stationary  solution  to  Eq.  (Bl)  as 

.  (B3) 

6u  ss  ^  Su^expdinBut)  . 

n  * 

This  assumption  leads  to  the  following  recurrence  relation 

(1  +  iin8cjT2)8un  =  i8Q"T2  (w(n  iy3  exp(i0')- w(o+lv3  exp(-i0") )  .  (B4) 

The  gain  of  the  in-quadrature  subharmonic  probe  field  is  proportional  to 
8ujexp(-i0").  From  Eq.  (B4)  we  find 

8uj  exp(-i0")  =  i8fTT2  w0  /  (1  +  *i8coT2)  .  (B5) 

The  right  hand  side  of  Eq.  (B5)  is  independent  of  0";  this  implies  the  in-quadrature 
subharmonic  probe-field  gain  is  independent  of  the  phase  difference  between  the 
oscillations  of  the  in-quadrature  probe  field  and  the  oscillations  of  the  strong  field. 
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Figure  Captions 

Figure  1.  The  applied  field  spectrum.  We  plot  the  spectrum  of  the  field  used  for  the 
calculation  as  vertical  arrows.  The  field  is  assumed  to  be  symmetrically 
placed  about  the  atomic  resonance,  <aatom.  The  atomic  resonance  is 
graphically  displayed  as  a  Lorentzian  lineshape.  The  field  isroade  up  of 
frequency  components  at  to  -  38u,  co  -  800,  gj  +  8u,  and  <o  +  38o.  The 
components  at  oa  4-  38co  and  o  -  38co  are  referred  to  as  the  strong  field  in 
the  text. 

Figure  2.  In-phase  subharmonic  probe-field  gain.  We  plot  the  in-phase 
subharmonic  probe-field  gain  as  a  function  of  8co,  for  various  values  of  the 
strong-field  intensity.  A  dip  in  the  center  of  the  gain  curve  develops  at 
low  intensity.  This  means  that  the  maximum  gain  occurs  for  nonzero 
modulation  frequencies.  The  time-averaged  strong-field  intensity,  Ij  = 
2(k£1)2T1  T2  =  0,  i,  11,  2, 2f,  31,  4,  41,  51,  6.  The  inversion  decay  time  Tj 

=  10  t2. 

Figure  3.  In-quadrature  subharmonic  probe-field  gain.  We  plot  the  in-quadrature 
subharmonic  probe-field  gain  as  a  function  of  8co,  for  various  values  of  the 
strong-field  intensity.  In  contrast  to  the  in-phase  subharmonic  probe- 
field  gain,  the  in-quadrature  subharmonic  probe-field  gain  develops  an 
antidip  near  the  center  of  the  gain  curve.  The  intensities  and  decay  times 
are  the  same  as  in  Fig.  2. 

Figure  4.  In-phase  subharmonic  probe-field  gain.  We  plot  the  in-phase 
subharmonic  probe-field  gain  as  a  function  of  8co,  for  larger  values  of  the 
strong-field  intensity.  Structure  develops  within  the  dip  that  formed  at 
low  intensity.  The  additional  structure  appears  at  subharmonics  of  the 
main  resonance.  The  time-averaged  strong-field  intensity,  It  = 
2(k£j)2T1T2  =  10,  20,  30,  40,  50,  60,  70,  80,  90,  100.  The  inversion  decay 
time  Tx  =  10  T2. 
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Figure  5.  In-quadrature  subharmonic  probe-field  gain.  We  plot  the  in-quadrature 
subharmonic  probe-field  gain  as  a  function  of  8co,  for  larger  values  of  the 
strong-field  intensity.  Additional  structure  appears  outside  the  antidip 
formed  at  low  intensity.  These  peaks  also  appear  at  subharmonics  of  the 
main  resonance.  The  time-averaged  strong-fieid  intensity,  Ij  = 
2(k£1)2T1T2  =  1C,  ^0,  30,  40,  50,  60,  70,  80,  90, 100.  The  inversion  decay 
time  Tx  =  10  T2. 

Figure  6.  Peak  positions  of  the  in-phase  subharmonic  probe-field  gain.  We  plot  the 
peak  positions  of  the  in-phase  subharmonic  probe-field  gain  as  a  function 
of  Rabi  frequency,  k^,  and  modulation  frequency,  8ui.  The  position  of  the 
resonances  in-modulation  frequency  varies  linearly  with  the  strong-field 
Rabi  frequency.  The  inversion  decay  time  Tx  =  10  T2. 

Figure  7.  Peak  postions  of  the  in-quadrature  subharmonic  probe-field  gain.  We  plot 
the  peak  positions  of  the  in-quadrature  subharmonic  probe-field  gain  as  a 
function  of  Rabi  frequency,  k^,  and  modulation  frequency,  8co.  The 
position  in  modulation  frequency  of  the  resonances  of  the  in-quadrature 
subharmonic  probe-field  gain  also  vary  linearly  with  strong-field  Rabi 
frequency.  The  slopes  of  the  lines  are  equal  to  the  slopes  of  the  lines  in 
Fig.  6.  The  decay  times  are  the  same  as  in  Fig.  6. 

Figure  8.  Strong-field  gain  curve.  We  plot  the  strong-field  gain  ns  a  function  of 
modulation  frequency,  8co.  The  maxima  occur  at  the  zeroes  of  the  zeroth- 
order  Bessel  function,  J0(2k£,  /  8co).  The  Rabi  frequency  of  the  strong  field 
is  =  5/T2.  The  inversion  decay  time  Tx  =  10  T2. 

Figure  9.  The  gain  curve  of  the  three  pairs  of  fields.  We  plot  the  (a)  strong-field 
gain,  (b)  in-phase  subharmonic  probe-field  gain  and  (c)  in-quadrature 
subharmonic  probe-field  gain  as  a  function  of  modulation  frequency,  8co 
for  the  same  parameters  as  Fig.  8.  The  absolute  maxima  of  the  strong- 
field  gain  curve  exceeds  the  subharmonic  probe  field  gains.  At  that 
modulation  fequency  the  strong-field  will  be  stable  to  the  growln  of 
subharmonic  probe  fields. 
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Figure  10.  The  influence  of  damping  on  the  in-phase  subharmonic  probe-field  gain. 

We  plot  the  in-phase  subharmonic  probe-field  gain  as  a  function  of 
modulation  frequency  for  three  different  ratios  of  T2  to  Tr  The  Rabi 
frequency  of  the  strong  field,  =  3/T2.  The  inversion  decay  times  for 
the  three  curves  are  (a)  Tx  =  £T2 ,  (b)  Tt  =  T2 ,  and  (c)  Tt  =  lOf^. 

Figure  11.  The  influence  of  damping  on  the  in-quadrature  subharmonic  probe-field 
gain.  We  plot  the  in-quadrature  subharmonic  probe-field  gain  as  a 
function  of  modulation  frequency  for  three  different  ratios  of  Tt  to  T2.  The 
Rabi  frequency  of  the  strong  field,  =  3/T2.  The  inversion  decay  times 
for  the  three  curves  are  (a)  Tj  =  £T2  ,  (b)  Tx  =  T2  ,  and  (c)  Tj  =  10T2. 
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NEEDS  FOR  BROADBAND  MATERIALS  AND  SYSTEMS 


•  Broadband  emitters  in  the  1  micron  band  for  either  flashlamp  or  diode- 
array  pumping. 

•  Broadband  absorber  in  the  800  nm  region  to  match  diode-array  emission, 
with  a  singie  or  broadband  emissions  in  the  1  micron  region. 

•  Non-linear  materials  to  provide  second  and  third  harmonic  generators, 
and  OPO  materials  for  the  visible,  the  near  I R,  the  mid  I R  and  the  far  I R,  at 
room  temperature. 

•  Broadband  emitters  (tunable  materials)  in  the  one-micron  band  based 
upon  the  Tunabe-Sugano  crystal  field  dependence  of  the  transition  metal 
energy  level  system. 

Example  of  Cr3  +  (a  d3  electron  system)  in  sapphire,  alexandrite. 

These  last  two  examples  are  not  suitable  because  they  lase  in  the  wrong 
band,  (around  750nm)  and  their  lasing  transition  lifetimes  are  too  short, 
(order  of  10  microsec  or  less). 

Ti  3  +  offers  a  possibility,  with  its  d'electronic  state,  to  provide  a  system 
with  little  or  no  deleterious  excited  staff  absorption. 

These  are  phonon  terminated  laser  transition  which  give  rise  to  the 
tunability  at  the  desired  wavelengths. 
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TANABE-SUGANO-TYPE  DIAGRAM  FOR  Ti3+ 
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LIST  OF  EQUIPMENT  AVAILABLE 

•  CW  dye  laser  system  including  a  20  W  Argon  laser  and  a  stabilized  ring  dye 
laser  (Coherent) 

•  Various  ancillary  measuring  equipment  such  as:  lock-in  amplifiers,  box-car 
integrators,  high-speed  storage  oscilloscope 

•  Various  computers  and  peripherals 

•  Optical  light  chopper 

•  X-ray  fluorescence  spectrometer  (Kevex) 

•  X-ray  generator  (Philips)  for  diffraction  analysis 

•  Visible  and  near  IR  interferometer  (Zygo) 

•  Pulsed  double-beam  spectrofluorometer  (Spex)  with  detectors  both  in  the 
visible  and  near  IR 

•  Streak  camera  and  monochromator  (Hamamatsu)  with  various  ancillary 
components 

•  Vacuum  pumping  system  (Baltzers) 

•  Single-crystal  x-ray  diffractometer  (Nicolet) 

•  Fume  hood 

•  Pulsed  Nd  dye-laser  system:  1  Joule  pump,  doubles  and  triples  converters 

•  Various  Newport  optical  tables  with  pneumatic  support 

•  Various  He-Ne  lasers 

•  Various  imaging  viewers  ramping  from  SI  photocathode  to  F.R  density 

•  Low-temperature  de _ (down  to  1Q0K 
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Transition  metal  ions  with  the  d1  or  d9  electronic 
configuration,  such  as  Ti3+,  V4+  and  Cu2+,  are  of  particular 
interest  for  tunable  solid  state  lasers.  In  these  ipns,  the 
existence  of  basically  only  one  high-lying  excited  state 
determines  the  absence  of  internal  transitions  that  could 
produce  excited  state  absorption.  Ti3+  in  corundum  is  an 
example  of  a  successful  system.  V4+  in  the  same  host  looks 
also  very  promising.1 

In  oxide  hosts,  vanadium  can  exist  in  a  number  of  ionic 
states  ranging  from  V2+  to  V5+.  V3+  is  by  far  the  most  stable 
ionic  state  of  vanadium  and,  because  of  this,  it  is  difficult 
to  obtain  significant  concentrations  of  V4+  in  this  host.  V4+ 
has  been  found  to  be  present  in  small  concentrations  in 
corundum  crystals  grown  by  the  Vemeuil  method.  Annealing 
under  oxygen  produces  the  conversion  of  some  V3+  into  V4+.2 
Gamma  radiation  is  also  effective  for  converting  V3+  into  V2+ 
and  V4+. 3  Codoping  with  magnesium  for  charge  compensation 
could  have  the  effect  of  increasing  even  further  the  ratio  of 
V4+  to  V3+  in  the  as  grown  crystals.1 


Crystals  of  vanadium-doped  corundum,  codoped  with 
magnesium,  were  grown  using  the  laser-heated  pedestal  growth 
(LHPG)  method.5  A  250  W  CO2  laser  allowed  us  to  grow  crystals 
up  to  3.5  mm  in  diameter.  The  LHPG  method  is  an  ideal  tool  for 
the  study  of  the  effects  of  growth  atmosphere  on  the  ionic 
valence  of  polyvalent  dopants,  such  as  vanadium,  since  it  is 
free  of  the  constraints  imposed  by  the  crucible  material  that 
limits  the  range  of  the  oxygen  fugacity  that  is  tolerable  in 
the  growth  atmosphere.  The  atmospheres  used  ranged  from  inert 
(low  oxygen  fugacity)  to  strongly  oxidizing  (about  47%  oxygen 
in  argon) .  An  atmosphere  with  18%  oxygen  gave  the  largest 
ratio  of  V4+  to  V3+  in  our  crystals,  as  judged  from  the 
absorption  spectrum,  when  the  feed  rods  were  prepared  starting 
with  an  equimolar  amount  of  magnesium  and  vanadium.  With 
higher  Mg/V  ratios,  the  optimum  was  found  to  be  at  lower 
oxygen  concentrations  (at  about  8%  of  oxygen  for  feed  rods 
prepared  with  a  starting  composition  containing  twice  as  many 
moles  of  magnesium  as  vanadium) .  The  final  composition  of  the 
crystals  could  not  be  determined  accurately  with  the 
analytical  method  available  in  our  laboratory  (x-ray 
fluorescence) .  The  vanadium  concentration  in  the  crystals  is 
certainly  much  lower  than  the  initial  concentration  in  the 
powder  mixture  as  a  significant  amount  of  vanadium  was  lost 
during  the  sintering  process  used  in  the  preparation  of  the 
feed  rods.  Further  losses  took  place  during  the  growth  of  the 
crystals. 

The  absorption,  fluorescence  and  excitation  spectra  for 
sample  #50  (size:  1.4  mm  x  1.6  mm  x  7.6  mm)  are  shown  in 
Figure  1.  These  spectra  were  obtained  at  room  temperature  using 
unpolarized  light  and  unoriented  crystals.  Their  main  features 
coincide  with  those  reported  in  the  literature. 1-4  The  orange 
emission  of  V4+  in  corundum  extends  from  520  nm  to  800  nm 
peaking  at  about  640  nm.  The  position  of  the  absorption  bands 
can  be  inferred  from  the  absorption  and  excitation  spectra. 
The  two  humps  peaking  at  425  nm  and  485  nm  can  be  identified 


as  related  to  the  V4+  intrinsic  absorption,  which  corresponds 
to  the  2T2g_2Eg  transition.  The  double  peak  characteristic  is 
due  to  the  dynamic  Jahn-Teller  effect  acting  on  the  2Eg  level 
(see  Ref.  4  and  references  therein).  The  features  with  maxima 
at  250  nm  and  320  nm  are  possibly  due  to  charge  transfer 
processes.  The  absorption  spectrum  of  sample  #37  (size:  2.0  mm 
x  2.5  mm  x  12.1  mm),  which  was  doped  only  with  vanadium  and 
grown  in  argon,  is  shown  for  comparison.  This  sample  contains 
mostly  V^+  which  has,  in  the  visible,  two  strong  absorption 
bands  with  maxima  at  400  nm  and  570  nm.  From  the  relative 
strength  of  the  absorption  bands  present  in  sample  #50  it  can 
be  concluded  that  most  of  the  vanadium  ions  in  this  sample  are 
in  the  4+  valence  state. 

The  fluorescence  lifetime  versus  temperature  was  also 
measured.  The  results  are  shown  in  Figure  1-d.  At  '  room 
temperature,  the  fluorescence  lifetime  is  1.2  j*s.  Thermal 
quenching  becomes  important  above  room  temperature.  From  these 
measurements,  and  the  width  and  position  of  the  emission 
spectrum,  the  resultant  value  for  the  stimulated  emission 
cross  section  is:  ase  *  5.2  10“19  cm2. 


References 


1.  V. A.  Sandulenko  and  E . A.  Sidorova,  Opt.  Spectrosc.  (USSR) 
6fi,  2  (1986). 

2.  J.Y.  Wong,  M.J.  Berggren  and  A.L.  Schawlow,  J.  Chem. 
Phys.,  ±2,  835  (1968). 

3.  M.  Villedieu,  N.  Devismes  and  A.M.  Goer,  J.  Phys.  Chem. 
Solids,  1&,  1063  (1977). 

4.  B.  Champagnon  and  E.  Duval,  J.  Phys  C:  Solid  State  Phys., 
12,  L425  (1979)  . 

5.  R.S.  Feigelson,  in  Tunable  Solid  State  Lasers,  ed  P. 
Hammerling,  A.B.  Bugdor  and  A.  Pinto  (Springer-Verlag, 
Berlin  Heidelberg,  1985),  p.  129. 


Optical  Density 


CECOM  CENTER  FOR  NIGHT  VISION  AND  ELECTRO-OPTICS 
CHARACTERIZATION  OF  THE  OPTICAL  SPECTRA  OF  EuA103:Ti  AND  GdA103:Ti 


1 


Characterization  of  the  Optical  Spectra  of  EuAlQ3:Ti  and 

GdA103 :Ti 


Larry  Merkle3,  Horacio  Verdun,  Uwe  Brauch*3, 
German  de  la  Fuente  and  Edward  Behrensc 
Fibertek,  Inc. 

510-A  Herndon  Parkway,  Herndon,  VA  22070 


Toomas  H.  Allik 

Science  Applications  International  Corporation 
1710  Goodridge  Drive,  McLean,  VA  22102 


Albert  Pinto  and  Wayne  Hovis 
U.S.  Army  Night  Vision  and  Electro-Optics  Center 
Fort  Belvoir ,  VA  22060 


aParmanent  address:  Department  of  Physics,  University  of 
Arkansas,  Fayettville  AR  72701 

^Present  address:  DFVLR  Institute  Fur  Technische  Physik, 
Pfaf fenwaldring  38-40,  D-7000  Stuttgard  80  Federal  Republic 
of  Germany 

cPresent  address:  Department  of  Physics,  Oklahoma  State 
University,  Stillwater  OK  74078 


2 


Abstract 


Two  oxide  perovskites,  EuAl03  and  GdA103 ,  doped  with 
titanium,  have  been  grown  by  the  laser-heated  pedestal 
growth  method.  The  crystal  structure  of  EuAlO-j  is  found  to 
be  different  than  that  reported  by  other  workers.  The 
absorption  spectra  of  EuAlC3:Ti  and  the  lack  of  fluorescence 
from  this  material  indicate  that  the  titanium  enters  the 
lattice  as  Ti4+.  On  the  other  hand,  the  spectra  of 
GdA103:Ti  are  attributable  to  Ti3+.  The  difference  between 
EuA103  and  GdAl03  can  be  attributed  to  the  difference  in 
third  ionization  potentials  between  europium  and  gadolinium. 
The  spectra  of  GdA103 : Ti  indicate  that  it  has  potential  as  a 
tunable  laser  material. 
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I.  Introduction 

The  trivalent  titanium  ion  placed  in  an  octahedral 
field  is  very  promising  for  tunable  solid  state  laser 
applications  due  to  its  strong,  broad  2eg-2t2g  emission 
band.  The  fluorescence  can  lie  in  the  long  wavelength 
visible  range  or  the  near  infrared,  depending  on  the  crystal 
field.  Its  [Ar]3d^  configuration  permits  no 
intraconf igurational  excited  state  absorption  transitions, 
which  is  a  significant  problem  in  some  other  solid  state 
systems.  1,2  Sapphire  doped  with  Ti3'4'  has  proven  to  be  a 
successful  tunable  laser  material.3'4 

There  is  significant  interest  to  investigate  Ti3+  in 
other  hosts,  in  which  laser  action  may  occur  in  spectral 
regions  other  than  that  of  Ti3  +  in  sapphire,  and  in  which 
crystal  growth  and  mechanical  and  thermal  properties  may  be 
optimized  for  laser  applications.  The  oxide  perovskites, 
wi^h  general  chemical  formula  AB03,  contain  a  large  A  ion 
occupying  a  twelvefold  coordination  site  and  a  smaller  B  ion 
in  the  sixfold  (octahedral)  site.  In  these  materials,  small 
transition  metals  tend  to  enter  the  B  site.5  Perovskites 
form  a  promising  class  of  laser  hosts  since  no  tetrahedrally 
coordinated  cation  sites  exist  whose  lower  crystal  field 
could  give  unwanted  absorption  at  the  sixfold  coordinated 
ions'  emission  wavelength. 
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The  work  reported  here  is  pair  of  a  systematic  study  of 
the  properties  of  Ti3+  in  perovskite  and  other  hosts  grown 
by  the  laser-heated  pedestal  (LHP}  growth  method.6  This 
method  permits  the  rapid  growth  of  single  crystals  of  a  wide 
range  of  compounds  so  that  more  extensive  studies  of  host 
dependent  optical  properties  can  be  undertaken  than  by 
Czochralski  or  Bridgman  techniques.  In  this  paper,  the 
optical  properties  of  EuAl03:Ti  and  GdA103 : Ti  are  reported. 
Section  II  reviews  the  pedestal  growth  technique  and  the 
crystal  structure  of  EuAlG-j  ,  and  Section  III  presents  the 
optical  properties  of  both  materials.  Interpretation  of  the 
data  is  discussed  in  Section  IV,  with  concluding  remarks 
following  in  Section  V. 

II.  Crystal  Growth  and  Properties 

The  crystals  for  this  study  were  grown  by  the  LHP 
growth  method  using  a  250-watt,  C0-.  laser.  Sintered  ceramic 
compacts  were  obtained  from  starting  stoichiometric  amounts 
of  A1203 (99.999%,Cerac) ,  Gd203 ( 99 . 999% , Cerac) , 

Eu203  (99 . 995% , Cerac) ,  and  Ti203  (99 . 8% , Cerac)  powders.  The 
starting  powdered  components  were  milled  in  acetone  for  ca . 
30  minutes  inside  a  sintered  corundum  jar  ball  mill 
(Brinkman  mixer  mill  model  MM2),  and  subsequently  dried  in 
an  oven.  The  dried  powder  was  then  placed  into  a  stainless 
steel  die  assembly  (28  mm  diameter)  and  a  compact  was  formed 
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at  a  force  of  14  tons  during  a  period  of  30  minutes.  The 
powder  compact  was  then  sintered  for  12  hours  at  1100°C 
(EuA103)  or  at  1500  3C  (GdA103)  in  air.  The  hard,  sintered 
compact  was  cut  into  rectangular  feed  nds  (of  size  ca .  2.5 
x  2.5  x  26  mm3)  and  mounted  cn  a  graphite  holder  fitted  into 
the  bottom  chuck  of  the  LHP  growth  apparatus.  A  platinum 
wire  was  attached  to  the  top  chuck  of  the  apparatus  and  used 
as  a  primer  to  initiate  the  nucleation  process.  The  LHP 
apparatus  is  comprised  of  a  CW,  250-watt,  C02  laser  adapted 
to  a  growth  chamber  capable  of  handling  oxidizing  and 
reducing  environments.  A  top  metal  chuck  inside  the  chamber 
is  used  for  mounting  primers  or  seeds,  while  a  bottom 
counterpart  is  used  for  mounting  ceramic  feed-rod 
precursors.  The  laser  beam  is  focused  onto  the  feedrod, 
which  is  rotated  in  order  to  evenly  distribute  the  heat 
generated  by  the  laser  throughout  the  molten  zone.  The 
rotating  seed  or  primer  is  then  inserted  into  the  melt  and 
nucleation  is  allowed  to  take  place.  After  several  minutes, 
the  crystal  is  pulled  from  the  melt  by  lifting  the  upper 
chuck  while  feeding  at  a  slower  rate  with  the  lower  chuck. 
The  resulting  polished  crystals  are  typically  1  to  2  mm  in 
diameter  and  20  to  50  mm  long.  This  growth  process  has 
also  been  recently  described  as  the  laser  floating  zone 
melting  method,7  and  has  been  adapted  for  the  growth  of 
materials  with  different  applications.8  The  specific  growth 
parameters  reported  here  are  summarized  in  Table  I.  Because 
both  titanium  and  europium  are  heterovalent ,  the  influence 
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of  the  growth  atmosphere  on  the  valence  state  must  be 
considered.  For  this  reason,  some  samples  were  grown  in 
air,  and  others  in  99.95%  pure  nitrogen. 

Single  crystal  x-ray  diffraction,  performed  on  a 
Nicolet  R3m/u  diffractometer,  was  employed  to  characterize 
the  EuAlO-j  samples.  Evaluation  of  several  axial  photographs 
reviel  that  structure  has  tetragonal  P  symmetry,  rather  than 
the  orthorhombic  structure  reported  by  others.9  The  unit 
cell  parameters  of  this  crystal  at  room  temperature  are 
aQ  =  0.3736  nm  and  cQ  =  0.7466  nm.  Since  the  pedestal 
growth  method  creates  large  thermal  gradients  along  the 
sample,  and  perovskites  have  the  propensity  for  structure 
transformations  in  narrow  temperature  ranges,10  it  is 
possible  that  the  observed  structure  represents  a  high 
temperature  phase  of  EuA103  "frozen  in"  at  room  temperature. 
This  structure  is  a  distortion  of  the  cubic  perovskite 
structure  with  a  cubic  pseudocell  size  equal  to  0.374  nm. 

III.  Optical  Spectroscopy 

Optical  spectra  were  taken  using  a  Spex  F222  and 
Perkin-Elmer  Lambda  9  and  MPF-66  spectrometers.  These 
instruments  permitted  the  measurement  of  absorption, 
fluorescence  and  excitation  spectra.  The  absorption 
spectrometer  is  equipped  with  a  4X  beam  condenser  (Harrick 
Scientific  Corporation,  Ossining,  NY).  Considerable  care 
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was  taken  to  discriminate  against  spurious  signals  such  as 
filter  fluorescence.  Fluorescence  lifetime  data  were  taken 
using  pulsed  excitation  by  the  second  and  third  harmonics  of 
a  Quantel  Q-switched  Nd:YAG  laser.  Detection  was  by  a 
cooled  photomultiplier  coupled  to  a  quarter  meter 
monochromator  or  by  a  silicon  photodiode  with  color  filters. 
Low  temperature  experiments  were  carried  out  by  placing  the 
samples  on  the  cold  tip  of  a  CTI  Model  21  cryogenic 
refrigerator. 

The  absorption  spectra  of  two  EuAlO-j  samples  grown  in 
air  are  shown  in  Figure  1  Absorption  in  the  undoped 
samples  for  wavelengths  longer  than  300  ran  is  attributable 
to  Eu3+  transitions,  Fluorescence  from  undoped  samples  also 
exhibits  only  transitions  of  this  ion.  No  evidence  of  Eu2^ 
absorption  or  emission  has  been  found,  unlike  the  results 
reported  by  Jaffe.11  The  room  temperature  fluorescence 
lifetime  of  the  Eu3  +  5Dq  state  in  the  undoped  sample  grown 
in  air  is  2.5  microseconds  for  355  nm  excitation  and  4.5 
microseconds  for  532  nm  excitation.  For  an  undoped  sample 
grown  in  nitrogen,  the  corresponding  values  are  6.5  and  17 
microseconds,  respectively.  In  this  sample,  lifetimes  were 
also  measured  at  low  temperature,  about  20  K,  and  were  found 
to  be  unchanged  from  the  room  temperature  values.  All  these 
lifetimes  are  far  shorter  than  observed  for  systems  with 
dilute  concentrations  of  Eu3+,  which  are  about  2  to  3 
milliseconds  even  in  glass.12  This  must  indicate  that  in 
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EuA103  the  lifetime  is  limited  by  interionic  processes  such 
as  energy  migration  to  traps.  The  difference  in  lifetimes 
between  the  two  samples  suggests  that  such  traps  are  more 
common  in  the  sample  grown  in  oxidizing  atmosphere. 

The  Eu3+  fluorescence  has  also  been  observed  in  lightly 
doped  EuAl03:Ti  samples.  Its  lifetime  at  room  temperature 
under  355  nm  excitation  is  0.5  microseconds  which  is  much 
shorter  than  in  the  undoped  samples.  This  suggests  that 
energy  transfer  from  the  Eu3  +  is  increased  in  the  doped 
samples . 

As  Figure  1  shows,  the  introduction  of  titanium  into 
EuA103  changes  its  absorption  spectrum  dramatically.  The 
strong,  very  broad  absorption  band  in  the  visible  and  near 
infrared  is  not  consistent  with  that  expected  for  Ti3  +  in  an 
octahedral  field.  In  addition,  doping  introduces  a  strong 
ultraviolet  band  between  about  400  nm  and  the  absorption 
edge  of  undoped  EuA103 .  The  strengths  of  both  these  bands 
increase  in  rough  proportion  to  the  Ti  concentration  in  the 
starting  materials,  as  shown  in  Figure  2.  Since  the  samples 
used  are  of  different  thickness,  the  absorbance  values  at 
345  nm  (indicative  of  the  strength  of  the  ultraviolet  band) 
and  at  650  nm  (the  visible  band)  are  normalized  by  the 
absorbance  of  the  Eu3+  peak  near  466  nm  in  each  sample. 

Note  that  two  samples  were  grown  with  0.48%  (atomic)  of  Ti 
in  the  starting  material,  and  that  both  absorption  bands  are 
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stronger  for  the  sample  grown  in  N2  than  for  that  grown  in 
air. 


No  fluorescence  has  been  observed  from  the  doped 
samples  except  for  weak  Eu3  +  emission  in  the  most  lightly 
doped  samples.  Taken  together  with  the  absorption  spectra, 
it  is  apparent  that  little  if  any  of  the  titanium  has 
remained  in  the  trivalent  state  in  these  samples.  The 
visible  band  does  not  correspond  to  intra-ionic  transitions 
in  other  charge  states  of  Ti,  but  its  strength  may  be 
indicative  of  a  charge  transfer  transition. 

Samples  of  GdAl03:Ti3+  have  also  been  grown  and 
sttidied.  In  these  samples  the  absorption  and  fluorescence 
in  the  visible  region  is  consistent  with  the  presence  of 
TiJ  .  The  excitation  and  fluorescence  spectra  are  presented 
in  Figure  3  for  a  sample  grown  in  nitrogen  with 
0.48  weight  %  Ti  added  to  the  melt.  The  pair  of  absorption 
peaks  at  445  and  490  nm  may  be  attributed  to  transitions 
from  the  2t2g  to  the  2eg  state,  using  the  familiar  labels 
which  would  apply  in  a  truly  octahedral  environment,  with 
the  latter  level  split  by  the  Jahn-rTeller  effect  as  well  as 
by  non-octahedral  components  of  the  crystal  field. 

Similarly  the  fluorescence  band  position  and  its  room 
temperature  lifetime  of  about  5  microseconds  are  consistent 
for  a  transition  from  the  lower  2eg  state  to  the  2t2g 
At  a  temperature  of  12  K  the  fluorescence  lifetime  was 


state . 


measured  to  be  14  microseconds,  thus  implying  a  quantum 
efficiency  of  36%  at  room  temperature.  The  positions  of 
these  transitions  indicate  that  the  crystal  field  at  the 
impurity  site  is  approximately  20%  stronger  in  GdA103  than 
in  A1203.3  The  change  expected  in  the  crystal  field 
approximation  can  be  estimated  by  noting  that  in  Al2c>3  each 
A1  has  three  oxygen  neighbors  at  a  distance  of  0.166  run  and 
three  at  0.197  run.  Whereas  in  GdA103,  assuming  for 
simplicity  the  cubic  perovskite  structure,  each  A1  has  six 
oxygen  neighbors  at  a  distance  of  0.1855  nm.13  The  familar 
inverse  fifth  power  distance  dependence  of  the  crystal  field 
thus  predicts  a  16%  stronger  field  in  GdA103  than  in  A1203, 
in  reasonable  agreement  with  experiment.  The  excitation 
peak  at  355  nm  does  not  correspond  to  any  transition  within 
the  3d1  configuration  of  Ti3  +  .  It  may  be  due  to  absorption 
by  another  impurity  which  then  transfers  its  energy  to  Ti, 
or  in  view  of  the  discussion  of  section  IV,  perhaps  to 
charge  transfer  from  Ti3+  to  Gd3  +  .  In  any  case,  it  is  clear 
that  the  titanium  dopant  enters  the  GdA103  lattice  in  the 
trivalent  state,  in  contrast  to  the  situation  for  EuA103. 

IV.  Discussion 

The  striking  difference  between  the  behavior  of 
titanium  in  EuA103  and  GdA103  is  probably  not  attributable 
to  crystal  field  effects,  as  the  ionic  radii  of  Eu3+  and 
Gd3+  are  quite  similar,  0.102  -  0.106  nm  and  0.100  - 


0.104 
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nm  respectively.14  Rather,  it  is  of  interest  to  note  that 
the  third  ionization  potentials  of  these  two  rare  earths  are 
quite  different,  24.92  eV  for  europium  and  20.63  eV  for 
gadolinium. 15  In  addition,  for  the  crystal  EuTi03,  the 
europium  ions  have  been  found  to  be  divalent  and  the 
titanium  ions  tetravalent.16,17  This  suggests  that  in 
EuA103:Ti,  the  dopant  ions  become  tetravalent,  each  making  a 
nearby  europium  divalent  to  conserve  charge.  The  ability  of 
Ti  to  remain  trivalent  in  GdA103  would  presumably  be  due  to 
the  smaller  third  ionization  potential  of  Gd,  inhibiting  it 
from  accepting  an  electron  from  Ti. 

A  simple  calculation  based  on  a  lattice  of  point  ions 
may  be  undertaken  to  test  the  suggestion  that  Ti  in  EuA103 
is  tetravalent  and  a  nearby  Eu  divalent.  Following  the  work 
of  Pedrini  et  al.18,  the  energy  difference  between  such  a 
state  and  a  situation  in  which  both  ions  are  trivalent  may 
be  evaluated  in  terms  of  the  ionization  potentials,  Madelung 
energies,  and  crystal  field  energies  of  the  two  ions  given 
in.Eq.  1. . 

E(Eu2++Ti4+)  -  E(Eu3++Ti3+)  =  I4 (Ti)  -  I3 (Eu)  -  EM(Ti  site) 

+  EM(Eu  site)  -  Ec (Ti-Eu)  -  Ecf(Ti)  (1) 

Here  I4(Ti)  is  the  fourth  ionization  potential  of  titanium, 
43.24  eV.19  I3(Eu)  is  the  third  ionization  potential  of 
europium,  EM  is  the  magnitude  of  the  Madelung  energy  at  the 


indicated  site,  Ec(Ti-Eu)  is  the  magnitude  of  the  Coulomb 
energy  of  attraction  between  the  extra  positive  charge  on 
Ti4+  and  the  extra  electron  on  Eu2+,  and  Ecf(Ti)  is  the 
crystal  field  contribution  to  the  energy  of  the  ground  state 
of  Ti3  +  .  The  ground  state  of  Ti4  +  is  nondegenerate  and  each 
valence  state  of  Eu  exhibits  only  a  small  crystal  field 
splitting,  so  that  the  crystal  field  contributions  for  these 
ions  are  neglected.  The  Madelung  energy  was  calculated  by 
standard  methods.20  To  simplify  the  calculation  of  the  EM 
and  Ec  terms,  the  lattice  has  been  approximated  by  that  of  a 
cubic  perovskite  with  lattice  constant  0.374  nm.  It  is 
assumed  that  the  Eu  which  accepts  an  electron  from  the  Ti  is 
one  of  the  nearest-neighbor  Eu  ions.  With  these 
assumptions,  EM(Ti  site)  is  45.95  eV,  E^ (Eu  site)  is 
22.57  eV  and  Ec(Ti-Eu)  is  4.46  eV.  The  crystal  field 
contribution  to  the  Ti3  +  ground  state  energy  is  estimated  by 
assuming  an  octahedral  environment,  so  that  ECf  =  -4Dq,  and 
was  determined  from  the  Ti3+  spectra  in  GdAl03  to  be  1.07 
eV.  In  this  way  it  is  estimated  that  a  Ti4+-Eu2+  pair  in 
Eu£103:Ti  is  lower  in  energy  than  a  Ti3+-Eu3+  pair  by  8.4 
eV.  Despite  the  crudeness  of  the  calculation,  this  result 
strongly  supports  the  attribution  of  the  lack  of  Ti3+  in 
EuA103:Ti  to  instability  of  the  trivalent  ion. 

Comparison  of  the  spectra  of  Ti  in  the  two  hosts 
suggests  that  the  above  calculation  has  overestimated  the 
energy  difference  between  the  two  charge  states.  Since  the 
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major  change  in  the  above  calculation  upon  switching  to  the 
GdAlO-j  host  is  the  4.29  eV  difference  in  1 3  of  the  rare 
earth,  this  calculation  predicts  that  Ti1  +  would  also  be 
unstable  in  this  host.  Thus,  the  presence  of  Ti3+  in  GdAl03 
may  suggest  that  the  energy  difference  between  the  two 
charge  states  has  been  overestimated,  and  that  the 
Eu3+-Ti3+  state  in  EuAl03:Ti  is  higher  in  energy  than  the 
Eu2+-Ti4+  state  by  no  more  than  about  4  eV. 

This  overestimation  is  surprising  since  a  similar 
approach  to  the  calculation  of  photoionization  energies  of 
divalent  rare  earth  dopants  has  generally  been  accurate  to 
within  an  electron  volt.21  The  larger  error  in  the  present 
case  makes  the  attribution  of  the  optical  properties  of 
EuA103:Ti  to  Ti4+  less  certain.  However,  the  error  may 
plausibly  be  attributed  to  the  crudeness  of  the 
calculations.  Not  only  has  an  idealized  (cubic)  version  of 
the  crystal  structure  been  used  to  calculate  the  Madelung 
energies,  but,  more  importantly,  no  attempt  has  been  made  to 
correct  these  energies  for  the  local  distortion  of  the 
lattice  around  the  impurity  site.  Although  no  data  are 
available  to  determine  this  distortion,  it  can  make  a 
sizable  contribution  to  the  Madelung  energies. 

The  above  discussion  suggests  plausible  interpretations 
for  the  strong  visible  and  near  infrared  absorption  band  and 
for  the  absorption  edge  near  400  nm  in  EuA103. 


The  visible- 
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infrared  band  may  be  attributed  to  a  charge  transfer 
transition  in  which  an  electron  is  moved  from  a  Eu2  +  ion  to 
the  nearby  Ti4+.  If  so,  the  energy  of  Eu2+-Ti4+  pair  must 
be  lower  than  that  of  a  Eu3+-Ti3+  pair  by  no  more  than  about 
one  eV.  The  400  nm  absorption  edge  lies  near  the  onset  of 
the  lower  absorption  band  of  Eu2  +  in  EuA103  reported  by 
Jaffe11  and  is  also  near  the  absorption  edge  of  Ti02,  in 
which  each  titanium  ion  is  coordinated  by  six  oxygens  much 
as  in  these  perovskites.  Thus,  the  ultraviolet  absorption 
may  be  due  to  Eu2+,  Ti4  +  ,  or  both.  It  should  be  noted  that 
no  evidence  of  Eu  fluorescence  has  been  observed  upon 
ultraviolet  excitation  of  EuA103:Ti.  However,  the  strong 
visible  absorption  band  overlaps  the  wavelength  range  in 
which  Eu  emission  is  expected,  so  that  efficient  energy 
transfer,  either  radiative  or  nonradiative ,  may  explain  this 
lack  of  fluorescence. 

The  excitation  and  fluorescence  spectra  of  GdA103:Ti 
permit  a  preliminary  consideration  of  its  potential  as  a 
laser  material.  The  five  microsecond  lifetime,  9  x  1013  Hz 
fluorescence  bandwidth  and  36%  quantum  efficiency  imply  a 
peak  stimulated  emission  cross  section  of  4  x  10~20  cm2. 

This  is  smaller  than  the  peak  cross  section  in  Al203:Ti  and 
comparable  to  that  of  typical  Cr3+-based  laser  materials. 

It  will  be  important  to  determine  the  nature  of  the  observed 
ultraviolet  excitation  band.  If  it  is  due  to  Ti3+  itself 
rather  than  to  some  inadvertant  impurity,  it  may  give  rise 
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to  excited  state  absorption.  The  wavelength  limit  of  the 
present  experiments  leaves  open  the  question  of  whether  such 
excited  state  absorption  would  interfere  with  laser 
emission . 

V.  Conclusions 

A  number  of  properties  are  important  in  determining  a 
material's  potential  as  a  solid  state  laser.  In  addition  to 
optical  properties  such  as  absorption,  excited  state 
absorption  and  stimulated  emission  cross  sections,  factors 
affecting  growth  such  as  size  and  charge  matching  must  also 
be  considered.  The  investigation  of  EuA103:Ti  presented 
here  demonstrates  the  influence  which  the  host  can  have  on 
yet  another  factor,  the  dopant's  valence  state.  It  is 
anticipated  that  similar  difficulties  would  arise  for  Ti 
doping  of  hosts  containing  Yb3+,  as  it  also  has  a  large 
third  ionization  potential.  Among  the  rare  earths,  only 
gadolinium,  and  perhaps  cerium  and  terbium,  has  a 
sufficiently  small  third  ionization  potential  and  wide 
enough  absorption-free  regions  in  the  visible  range  to  be  of 
interest  for  titanium-doped  laser  materials.  As  noted  in 
the  previous  section,  GdAlO-jiTi  does  indeed  appear  to  be  an 
interesting  system  for  further  study. 

The  point-ion  based  calculation  of  the  energy 
difference  between  the  two  possible  charge  states  of  the 
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dopant  and  a  nearest  neighbor  rare  earth  was  only  partially 
successful,  having  predicted  that  Ti4+  would  be  more  stable 
than  Ti3+  in  GdAl03  as  well  as  in  EuA103 .  Without  direct 
data  on  the  positions  of  the  ions  neighboring  the  impurity 
site,  it  is  difficult  to  know  whether  this  difficulty  is 
fully  attributed  to  the  crudeness  of  the  ionization 
calculations.  Based  on  the  known  charge  states  of  the  Ti 
and  Eu  ions  in  EuTi03,  however,  the  interpretation  of  the 
spectra  given  in  Section  IV  remains  the  most  plausible. 

Acknowledgments:  Both  the  Fibertek  and  SAIC  work  was 

sponsored  by  the  Center  for  Night  Vision  and  Electro-Optics. 


17 


References 


1.  P.  F.  Moulton,  Laser  Focus  19(5),  83  (1983). 

2.  M.  L.  Shand  and  H.  P.  Jenssen,  IEEE  J.  Quantum 
Electron.  OE-19 ( 3) .  480  (1983). 

3.  G.  F.  Albrecht,  J.  M.  Eggleston,  and  J.  J.  Ewing,  Opt. 
Commun.  52.(6),  401  (19  85). 

4.  Several  papers  on  this  system  are  found  in  Tunable 
Solid  -State  Lasers  II.  ed.  A.  B.  Budgor,  L.  Esterowitz,  and 
L.  G.  DeShazer,  (Springer-Verlag,  Berlin  Heidelberg,  1986). 

5.  J.  P.  van  der  Ziel  and  L.  G.  Van  Uitert,  Phys.  Rev.  B 
8(5),  1889  (1973). 

6.  R.  S.  Feigelson,  in  Tunable  Solid  State  Lasers,  ed.  P. 
Hammerling,  A.  B.  Budgor,  and  A.  Pinto  (Springer-Verlag, 
Berlin  Heidelberg,  1985),  p.  129. 

7.  K.  Dembinski,  J.  M.  Bassat,  J.  P.  Coutures  and  P. 

Odier,  J.  Mater.  Science  Lett.  6,  1365  (1987). 

8.  X.  P.  Jiang,  G.  W.  Quao,  J.  G.  Huang,  Y.  Yu,  M.  Jiang, 
Y.  L.  Ge,  Z.  Q.  Hu,  C.  X.  Shi,  Y.  H.  Zheo,  Y.  J.  Wang,  G.  Z. 


Su  and  Y.  E.  Zhuo,  Extended  Abstracts,  High-Temperature 
Superconductors  II,  MRS  Spring  Meeting  Symposium,  Reno,  NV. 
1988,  D.  W.  Capone  II,  W.  H.  Butler,  B.  Batlogg  and  C.  W. 
Chu,  EDS.  Materials  Research  Society,  Pittsburgh,  PA,  1988. 

9.  S.  Geller  and  V.  B.  Bala,  Acta  Cryst.  9,  1019  (1956). 

10.  O.  Muller  and  R.  Roy,  The  Major  Ternary  Structural 
Families .  (Springer-Verlag ,  New  York,  1974)  pp. 175-229. 

11.  P.  M.  Jaffe,  J.  Electrochem.  Soc.  117 (7) .  918  (1970). 

12.  X.  Gang  and  R.  C.  Powell,  J.  Appl.  Phys.  52(4),  1299 
(1985)  . 

13.  R.  Wyckoff,  Crystal  Structures.  (Wiley  Interscience, 
New  York,  1964)  vol.  2,  pp.7,392. 

14.  S.  K.  Dickinson,  Air  Force  Cambridge  Research 
Laboratories  Physical  Sciences  Research  Paper  439  (1970). 

15.  W.  Martin,  L.  Hagan,  J.  Reader  and  J.  Sugar,  J.  Phys. 
Chem.  Phys.  Ref.  Data  2(3),  771  (1974). 


16.  J.  Brous,  I.  FanKuchen,  and  E.  Banks,  Acta  Cryst.  6,  67 


19 


17.  T.  R.  McGuire,  M.  W.  Shafer,  R.  J.  Joenk,  H.  A. 

Alperin,  and  S.  J.  Pickart,  J.  Appl.  Phys.  22(3),  931 
(1966)  . 

18.  C.  Pedrini,  D.  S.  McClure,  and  C.  H.  Anderson,  J.  Chem. 
Phys.  70(11),  4959  (1979). 

19.  Handbook  of  Physics  and  Chemistry,  52nd  ed. ,  ed.  Robert 
Weast,  (Chemical  Rubber,  Cleveland,  1971)  . 

20.  M.  P.  Tosi,  Solid  State  Physics,  vol.16,  editors  F. 
Seitz  and  D.  Turnbull  (Academic  Press,  New  York,  1964) 
p.  13  . 


21.  C.  Pedrini,  F.  Rogemund,  and  D.  S.  McClure,  J.  Appl. 
Phys .  59(4),  1196  (1986). 


20 


Figure  Captions 


1.  Room  temperature  optical  absorption  of  two  EuAlO^ 
samples  grown  in  air.  Lower  trace:  Undoped,  sample  length 
is  1.14  mm.  Upper  dashed  trace:  Doped  sample  with  0.47% 
(atomic)  Ti.  Sample  length  is  0.71  mm.  (The  concentration 
refers  to  that  in  the  starting  material.) 

2.  Relative  absorbance  at  345  nm  and  at  650  nm  versus  x  in 
EuA11_x03 :Tix.  Each  value  is  a  ratio  of  the  absorbance  at 
the  stated  wavelength  to  the  absorbance  of  the  Eu3+  line  at 
466  nm. 

Squares:  Absorbance  at  650  nm  for  samples  grown  in  air. 

Asterisk:  Absorbance  at  650  nm  for  sample  grown  in  N2. 
Diamonds:  Absorbance  at  345  nm  minus  that  at  440  nm, 
samples  grown  in  air. 

Circle:  Absorbance  at  345  nm  minus  that  at  440  nm,  sample 

grown  in  N2 . 

3.  Room  temperature  fluorescence  and  excitation  spectra  of 
GdA103:Ti.  The  fluorescence  spectrum  is  corrected  for 
instrument  response. 
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